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suitable  for  computer  application  is  studied  with  special  emphasis  placed  on 
the  chemical  interactions  of  the  fuel  with  the  re-entry  environment.  General 
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tions  are  written  in  a  form  that  includes  effects  from  transpiration  cooling, 
heats  of  iusion,  evaporation,  dissociation,  recombination,  chemical  reactions, 
and  mass  accretion  (or  loss).  All  pertinent  parameters  are  evaluated  insofar 
as  available  data  will  permit.  The  rate  of  hydrogen  transpiration,  the  heat 
energy  added  to  the  fuel  element  due  to  recombination  of  the  transpired 
hydrogen,  and  the  vapor  pressure  of  molten  fuel  material  could  not  be  determined 
with  reasonable  accuracy  from  the  available  data.  It  is  recommended  that  these 
three  parameters  be  determined  through  laboratory  experiments  and  used  in  the 
model  described  to  determine  whether  or  not  the  fuel  elements  will  burn  up  during 
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I.  INTRODUCTION 


Ona  of  the  problems  attendant  with  the  use  of  fuel -rod  “carrying  nuclear  reactors 
in  satellite  applications  is  the  possibility  that  these  fuel  rods  will  re-enter  the  earth's 
atmosphere  without  burning  up  entirely  and  reach  the  earth  wholly  or  partially  intact, 
thereby  producing  radiation  hazards  in  populated  areas.  The  purpose  of  the  present  pro¬ 
gram  is  to  establish  a  mathematical  model  whereby  various  fuel  element  re-entry  conditions 
can  be  simulated  on  a  computer  to  determine  if  such  fuel  elements  will  indeed  burn  up  at 
sufficiently  high  altitudes  to  render  earth  surface  level  hazards  negligible. 

Several  previous  studies  on  the  fuel  element  materials  have  been  made,  as  is  evi¬ 
dent  from  the  bibliography  in  this  report.  These  data,  together  with  samples  of  fuel  element 
material  that  had  been  subjected  to  simulated  re-entry  tests,  formed  the  basis  of  the  present 
study. 

The  problem  of  describing  the  re-entry  behavior  of  nuclear  fuel  material  may  be 
divided  into  three  general  areas: 

(1)  Establishing  the  system  of  re-entry  equations; 

(2)  Delineating  the  parameters  that  must  be  known  in  order  to  solve  this 
system  of  equations;  and 

(3)  Through  chemical  analysis  on  the  samples  and  mathematical  analysis  of 
the  available  data,  determining  the  values  of  these  parameters.  Those 
parameters  that  are  not  constant  throughout  the  re-entry  regime  must  be 
described  in  analytical  form  for  simultaneous  solution  with  the  re-entry 
equations. 

The  system  of  equations  has  been  established  for  the  limited  case  of  re-entry  of  cy¬ 
lindrical  fuel  rods,  although  the  analysis  may  be  applied  to  the  more  general  case  of  re¬ 
entry  of  material  of  arbitrary  shape.  The  assumptions  used  in  deriving  this  system  of  equations 


1 


are:  stationary  planet  atmosphere;  planar  trajectory;  negligible  effect  of  net  lifting  forces 
on  the  trajectory;  and  uniform  wall  temperature  due  to  the  rod  spinning  about  its  cylindri¬ 
cal  axis.  The  parameters  required  for  the  solution  of  this  set  of  simultaneous  differential 
equations  have  been  defined  and  their  values  determined  as  accurately  as  the  existing  data 
will  allow. 

The  primary  limitations  on  this  study  are  due  to  the  inadequacy  of  the  available  data. 
Some  of  those  inadequacies  are  discussed  at  the  end  of  the  report  and  possible  measures  to 
overcome  their  shortcomings  are  suggested. 
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II.  EQUATIONS  OF  MOTION  AND  RE-ENTRY 


In  this  section,  the  trajectory  and  angular  motion  equations  will  be  formulated  to¬ 
gether  with  the  effects  of  ablation  and  accretion  on  these  equations.  The  evaluation  of 
the  various  parameters  involved  in  these  equations  is  the  subject  of  the  following  section. 

A.  RE-ENTRY  EQUATIONS 

Fi  ve  assumptions  are  employed  in  deriving  the  equations  of  motion  of  a  re-entering 
vehicle.  The  first  two  of  these  are: 

Assumption  1 .  The  planet  atmosphere  is  stationary . 

Assumption  2.  The  re-entry  trajectory  is  planar. 

The  following  notation  is  employed: 

rQ  =  radius  of  earth.  (This  can  be  considered  as  the  distance  from  the  center 
of  gravity  to  the  surface  of  a  reference  sphere  having  the  average  radius 
of  the  earth.) 

y  =  altitude  above  earth's  surface  (can  be  considered  as  altitude  above  the 
surface  of  the  reference  sphere) . 

x  =  co-ordinate  along  the  surface  of  a  sphere  from  the  point  of  re-entry, 
i  .e„  (x,y)  =  (0,  yg)  at  the  point  of  re-entry. 

r  =  distance  from  center  of  gravity  of  vehicle  to  center  of  gravity  of  earth. 

rc  =  radius  of  curvature  of  trajectory. 

8  -  trajectory  angle  cbove  horizon  of  sphere  of  radius  (r©  +  y) . 

ds  =  distance  along  trajectory  (ds8  =  dy8  +  dx8). 

V  =  velocity  of  vehicle  relative  to  stationary  atmosphere, 
p  =  density  of  atmosphere. 

A  =  cross  sectional  area  of  vehicle  ( id),  l  =  length,  d  =  diameter  of  the 
vehicle. 
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W  =  weight  of  vehicle  (at  earth's  surface),  and 
g  =  acceleration  of  gravity  (at  earth's  surface),  so  that 
W/£j  =  mass  of  vehicle . 

l 

La  =  lift  of  vehicle  =  F  PVaACi_* 

Cl  =  lift  coefficient. 

D  =  drag  of  vehicle  =  j  P Vs A Cq. 

Cq  =  drag  coefficient. 


t  =  time  (token  as  zero  at  point  of  re-entry,  i  .o.  (x,  y,  t)  =  (0,  yjr,  0)  at 
point  of  re-enrr/).  E  is  a  subscript  denoting  conditions  at  the  point  of 
re-entry. 

Other  symbols  will  be  defined  as  they  are  brought  into  use. 


Note  .hot  d£|  =  .  tgn  0E  _ 


Now,  it  is  appropriate  to  derive  a  relationship  between  rQ  and  rc .  To  do  this,  refer 
to  Figure  1 .  From  the  geometry,  it  follows  that 


ds  =  rcdY,  ds  =  rQ  sec  0do,  and  -d  0  *  dY  -  da. 


where  d  Y  is  the  elementary  angle  subtended  by  the  elementary  arc-length,  ds,  along  the  ve¬ 
hicle  trajectory,  and  do  is  the  elementary  angle  si4>tended  by  ds,  as  viewed  from  a  point  at 
the  center  of  gravity  of  the  earth.  Therefore, 

d0  _  dY  d<7  _  J_  cos  8 
"  ds  3T  ’  dT  ~  rc  rQ 


Thus,  the  following  identities  hold: 


JL  =  cos  9 


de 

IT 


ss 


de 

IT 


for  small  0  ,  and 
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\ 

i 

i 

|  Figure  1 .  Diagram  showing  the  relationship  between  coordinates  and  paraaeters 

■  at  a  point  on  a  re-entry  trajectory. 

i 

\ 
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Hence , 


d 0  _  I  d0 

cIT  V  ar  • 


1  _  cos  8  1  d  9 

rc  =  ro  “  V  ‘ 

For  forces  normal  to  the  flight  trajectory.  Figure  2  shows  that 


L  J  a  r- 


Assumption  3.  Since  the  net  effect  of  lift  forces  on  the  fuel  element  in  question  (a 
cylinder  whose  length  is  approximately  ten  times  its  diameter)  only  changes  the  length  of  the 
trajectory  by  a  small  fraction,  they  will  be  assumed  to  be  negligible;  hence. 


-  W 


[?] 


a  WV3 

cos  0  =  - - 

g  rc 


W  [~cos  8  1  d  0  “I 

-  -  -  v2!.—  -  V  17  J 


or 


E!®  =  V  cos  8  f±_  . 

W  L  r_ 


9ro  1 

WJ 


0) 


For  forces  along  the  flight  path,  Figure  2  shows  that 

D  +  W[i>n..-^V]. 

Making  use  of  the  definition  of  the  drag,  D,  and  performing  the  indicated  differentiation, 
this  last  expression  becomes 

ipN^ACo+wl^r]  sin0  =  -y  17"  y  ^7  *  (2) 

By  definition  of  coordinates  and  0  ,  it  follows  that 

g-7  =  V  cos  9  (3) 


gj-  “  V  sin  0  . 


6 


(4) 


Figure  2.  Diagrams  showing  the  relationships  between  coordinates 
and  vectors. 
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Figure  2.  Diagrams  showing  the  relationships  between  coordinates 
and  vectors. 
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These  four  equations  (eqs.  1 ,  2,  3,  and  4)  specify  the  re-entry  motion  of  the  fuel  elements 
in  terms  of  the  initial  conditions  and 

d  6  dV  dx  dy 

Wm 3t'3T'3t  • 

The  quantities  V,  r,  W,  and  8  can  be  calculated  from  the  foregoing  equations  if  the 
quantities  p  ,  Cp,  A,  and  <M/dt  can  be  specified  independently. 

The  value  of  p  can  be  determined  either  from  atmospheric  tables  or  from  the  approxi¬ 
mation: 

— 8  y 

P  =  ap0e  y 

where 

apQ  =  sea  level  density  (slugs/ft3), 

a  =  0.715,  and 

B  =  1/(24,800  ft). 

The  values  of  A  and  CD  must  be  determined  from  the  angular  motion  of  the  rod,  since 
they  depend  on  the  orientation  of  the  fuel  element  for  their  specific  values  during  any  instant 
of  time  (t,  t  +  dt).  The  dynamics  of  a  re-entering  rod  are  fairly  well  understood.  (See,  for 
example:  Raymond  and  Garber,1  Davey  and  Grigsby,8  and  others.) 


B.  ANGULAR  MOTION  EQUATIONS 

The  general  angular  motion  (for  a  rod  tumbling  in  the  plane  of  the  trajectory)  should 
be  governed  by  an  equation  of  the  form9  * 


*  F(«  .  A,  d  ,  da)/  Voo),  (4a) 

•-&'  <*> 

‘Raymond  assumes  the  static  restoring  and  damping  forces  are  respectively  proportional  to 
sin  2 <t>  and  cos80  .  AVCO's  calculation3  indicates  that  this  is  not  for  from  correct. 


8 


where; 


F  = 

co  = 

u>  - 
<t>  = 

I  = 
l  = 
d  = 

A  = 

Poo= 

Performing  the  indicated  differentiation  in  eq.  (6a),  making  use  of  eq.  (6c),  and  sub¬ 
stituting  an  analytical  fit  to  AVCO's3  calculations  on  the  static  restoring  and  damping  coeffi¬ 
cients  for  F,  the  angular  equation  of  motion  becomes* 

2 

+W^]=  1.4PooV^)Adsin20[l+2.17sin320]  -  71  P^V^Ad*!  sin*| .  (7) 

9 

Assumption  4.  It  has  been  assumed  that  rg  is  constant  in  deriving  eq.  (7) .  When 
weight  losses  and/or  dimensional  changes  are  such  that  rg  cannot  be  considered  constant,  the 
term  (2u)W Tg/g)  d  rg/dt  must  be  added  to  the  left-hand  side  of  eq .  (7) . 

The  two  additional  parameters  duu/dt  and  d$/dt,  introduced  by  eq.  (7),  can  be  deter¬ 
mined  by  the  simultaneous  solution  of  equations  (6b)  and  (7)  with  equations  (1)  through  (5), 
provided  that  the  initial  conditions  on  0  and  m  are  specified.  Due  to  the  relatively  small 

*Note,  AVCO's  curve  (Fig.  B-l-8  in  ref.  3)  could  be  read  point  for  point  into  a  computer 
program  and  used  in  place  of  the  right-hand  side  of  eq.  (7).  The  overall  accuracy  of  this 
analytical  fit  to  AVCO's  curve  is  within  5%. 


is  a  function  to  be  defined  in  eq.  (8). 
subscript  denoting  free  stream  conditions, 
angular  velocity  of  rod . 

angle  of  cylindrical  axis  of  rod  with  respect  to  end-flow  position. 
(0  =  rr/2  at  cross-flow;  0  =  0  for  end-on  flow  condition.)  0  is 
termed  the  pitch  angle. 

—  r®  =  moment  of  inertia,  with  rg  =  radius  of  gyration. 
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length  of  rod  (feet), 
diameter  of  rod  (feet). 

16. 

free  stream  density  (slugs/ft3). 
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thermal  expansion  of  the  fuel  element  material  considered  in  this  report,  the  effect  of  this 
expansion  on  rg  can  be  neglected  for  practical  purposes.  It  should  be  emphasized  that 
eq.  (7)  is  valid  only  for  cylindrical  vehicles.  If  other  vehicle  shapes  are  to  be  used,  the 
right-hand  side  of  this  relation  has  to  be  evaluated  for  the  particular  shape  in  question. 

C.  ABLATION  EQUATIONS 

In  order  to  discuss  accretion  and  ablation,  it  is  necessary  to  consider  the  heat  transfer 
to  the  rod,  the  temperature  of  the  wall,  etc .  The  quantity  that  can  be  discussed  with  the 
greatest  degree  of  confidence  is  the  rate  of  heat  transfer  to  the  fuel  element. 

t 

Assumption  5.  It  is  assumed  that  the  fuel  element  will  "spin"  about  its  cylindrical 
axis  so  that  the  wall  can  be  considered  to  be  at  a  uniform  temperature.  Because  of  minor 
surface  and/or  internal  irregularities,  especially  where  tumbling  is  also  taking  place,  this 
should  be  an  accurate  assumption  up  to  the  time  large  pieces  of  the  rod  have  been  ablated 
away.  That  is,  fairly  large  asymmetries  would  be  required  to  assure  the  stability  of  a  cylin¬ 
der  during  re-entry. 

Assumption  6.  It  will  be  assumed,  for  the  moment,  that  the  fuel  element  remains  at 
a  uniform  temperature  throughout.  This  assumption  will  be  examined  in  the  next  section, 
and  a  method  for  compensating  for  the  fact  that  this  assumption  is  not  strictly  correct  will  be 
presented . 

With  these  assumptions,  the  instantaneous  heat  balance  equation  (in  the  absence  of 
transpiration  cooling  or  ablation)  requires  that 

=  [CH  fbVoo  -  hw>  +  H=h  “  €oTw]^  *&"  *  (8) 

where  p  ,  V,  W,  and  g  are  as  defined  previously  (ao  refers  to  free  stream  value). 
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Tw  =  temperature  of  wall  (°R). 

CH  =  average  heat  transfer  coefficient  (averaged  over  one  cycle  of 
"spin")  -  dimensionless. 

he  =  stagnation  enthalpy  at  outer  edge  of  boundary  layer  (Btu/slug). 

hw  =  enthalpy  of  air  at  the  wall  temperature  (Btu/slug). 

e  =  emittance  of  the  wall  material  -  dimensionless. 

a  =  Stefan-Boltzmann  constant  =  4.8  x  10“13  Btu/ft8 -sec °R* . 

Hch  =  rate  of  chemical  heat  input  to  rod  due  to  chemical  reactions 
(Btu/ft8-sec) . 

S  =  surface  area  of  rod  (ft3) 

Cp  -  average  specific  heat  capacity  of  rod  (Btu/slug  -  °R),  and  may 
be  a  function  of  its  temperature  history. 

The  heat  transfer  coefficient,  C^j/  which  depends  on  the  re-entry  conditions  and  geome¬ 
try  of  the  vehicle  in  question,  can  be  obtained  from  published  literature,  ostensibly  from  papers 
by:  R .  J.  Vidal  and  C  .  E .  Wittliff;4  and  J.  C  .  Fay  and  F .  R .  Riddell .5  The  value  of  Tw  is 
obtained  from  the  initial  conditions  and  eq.  (8),  solved  simultaneously  with  the  preceding 
seven  equations.  The  emittance,  6  ,  and  the  average  heat  capacity,  Cp,  must  be  known  for 
every  value  of  Tw.  The  quantity  Hch  must  be  calculated  from  the  environmental  conditions  and 
the  chemical  interaction  of  the  rod  material  with  the  environment. 

The  change  in  weight  rate,  dW/dt,  may  be  subdivided  into  three  categories:  the  decrease 
in  mass  due  to  outgassing  of  the  rod,  the  decrease  in  weight  due  to  fusion  (melting)  and  evapora¬ 
tion,  and  the  increase  in  weight  due  to  chemical  interactions. 

In  the  outgassing  category  (dehydriding),  first  assume  that  no  outgassing  takes  place  un¬ 
til  a  certain  wall  temperature,  TWg,  is  reached,  and  that  Tw  =  TWg  until  all  the  gas  has  escaped. 
Then  the  heat  balance  equation  for  transpiration  cooling  requires  that 
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(9) 


1  dW 
g  "cF 


[Fh  poo Yx>^e  ~  ^wg)  4  *^ch  ~  c p^4wgJ £ 


deG,  Twg  HdeG  +  ri(hfl  -  hwg)  -  Fc  HC/  deG 


where 


—  ^j~|  =  rate  of  change  in  mass  of  the  rod,  in  slugs/second, 

9  ldeG,TWg  due  to  outgassing . 

^wg  =  the  outgassing  temperature . 

HdeG  =  ener9y  required  to  release  one  slug  of  gas  (Btu/slug). 

^c,  deG  =  heat  energy  released  in  the  combustion  or  re-combination 
of  the  transpired  gas  (Btu/slug). 


r) 

Fc 


=  transpiration  coefficient  (dimensionless). 

=  fraction  of  the  combustion  energy,  Hc  deG,  that  is 
absorbed  by  tlte  fuel  rod  (dimensionless) 


hwg  =  enthalpy  of  air  at  the  temperature  Twg. 

Equation  (9)  neglects  the  rate  of  increase  (or  decrease)  of  weight  due  to  Hch.  This 
will  be  treated  separately.  The  last  term  in  the  denominator  of  eq.  (9)  accounts  for  the  effects 
due  to  possible  changes  in  the  state  of  an  absorbed  gas  after  it  has  been  released  from  fuel  rods 
such  as  the  combustion  or  re -combi  nation  of  the  hydrogen  as  it  is  released.  This  relation  is 
valid  so  long  as  the  denominator  is  positive. 

The  value  of  "H  has  been  studied  by  Rose  and  Allenhartz6  and  by  Georgiev  7  It  is  given 


by 


r  29  i 
1  -  0.60LMJ.GJ 


Ojas 


(10) 


where  ^JeG  '*  the  molecular  weight  of  the  outgassed  material .  This  will  be  discussed  in  more 
detail  in  the  following  pages.  In  case  more  than  one  molecular  weight  is  involved,  I* 

the  average  molecular  weight,  i  .e . ,  the  weight,  in  grams,  of  a  mole  of  the  transpiring  gas 
mixture 
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In  the  case  where  the  temperature  does  not  remain  constant  during  the  blowing  of  the 
outgassing  material,  the  situation  becomes  more  awkward  to  treat  mathematically.  It  can  be 
handled  in  a  formal  manner  when  rod  integrity  (shape)  is  maintained .  According  to  the  theory 
of  transpiration  cooling,*'7  the  heat  flux  through  the  wall,  q^,  with  transpiration  is  less  than 
the  value,  q0,  without  transpiration  cooling  by  the  factor  -n (he  -  hw)drr/dt.  The  quantity 
dm/dt  is  {ust  the  rate  of  mass  loss  per  unit  area;  hence. 


1  dW  I 

where  — -ir*  Is  the  magnitude  of  the  weight  loss  rate  per  unit  area  due  to  the  evolution 
l9  *deG 

of  the  transpiring  gas.  Now  the  heat  balance  condition  requires  that 

6  |1  dWj  WCp  dTw 

+  ftch  -  (H,  *  Hv)|-  -^|deG  ♦  1/  -af  +  £!,Ti- 

Here,  Hf  and  Hv  are  respectively  the  heats  of  fusion  and  vaporization  of  the  material  that  is 
producing  the  evolving  gas  (Btu/slug).  Eliminating  qtr  from  these  two  equations  and  noting 
that  qc  =  “  Kv)/  th®  expression  for  the  rate  of  change  in  the  total  vehicle  mass. 


due  to  transpiration,  is 


1  dwl 

‘•■*1 


[  <UVa>CH(h.  -  ^>+Hch-e0Ti.M]s 


Hv  +  «f  +  »l(h.  -  K) 


This  equation  is  to  be  used  in  place  of  eq.  (9)  when  Tw  is  not  constant  during  outgassing. 

There  may  be  an  accretion  of  mass  associated  with  the  chemical  reactions  whose  heat 
release  (or  absorption)  rate  is  Hch.  If  M,-  and  are  respectively  the  masses  of  the  rod  ma¬ 
terial  and  the  accreted  compound  in  slugs/gram-mole,  then,  clearly,* 


1  dwl  _  rHchl 

•  *JL'  Ll^j 


(Mr  -  Mch)$, 


*Hch  is  negative  for  exothermic  reactions,  but  |rlcfo  always  has  the  opposite  sign  of  H^. 


(12a) 


where  Hcj,  is  the  heat  of  the  reaction  M,.  +  O  +  N - ¥  Mch  in  Btu/gram-mole.  (The  O  and 

N  may  be  molecular,  partially  ionized  or  atomic  in  this  reaction.)  If  there  are  several  re¬ 
actions  taking  place,  a  summation  over  Hch,  Hc^,  and  Mch  mus*  be  employed.  When  eq.  (12a) 
is  added  to  eq .  (1 1),  the  net  rate  of  mass  loss  (or  gain)  is  obtained .  These  equations  are  valid 
so  long  as  the  accreted  products  remain  on  the  rod.  When  (or  if)  the  accreted  products  start  to 
flakrf  off  of  the  rod,  a  different  treatment  is  required. 

In  view  of  the  fact  that  both  ^1  and  appear  in  eq.  (1 1),  an  equation  similar 

dt  IdeG  dt 

to  eq.  (12a)  must  be  deduced  for  separating  these  two  parameters.  For  the  present,  this  independ¬ 
ent  relationship  will  be  written  as 

-  1  ^  =  G(T„)  (12b) 

9  dt  deG 

Finally,  the  case  of  fusion  and  evaporation  of  the  rod  will  be  considered.  This  process 
is  considered  to  take  place  at  the  fusion  temperature,  Twf,  and  a  certain  fraction,  Fv,  of  the 
material  is  assumed  to  evaporate  and  produce  transpiration  cooling.  By  similarity  with  eq .  (9), 
this  vehicle  mass  change  rate  is  obviously  given  by 


It  should  be  noted  that  eqs.  (1 1)  and  (12)  are  to  be  solved  simultaneously  with  eqs.(l) 
through  (8).  In  a  case  where  eqs.  (9)  and  (13)  are  applicable,  the  first  eight  equations  are 
followed  until  the  respective  wall  temperatures  reach  TWg  and  Twf,  and  the  temperature  is 
then  held  constant  (at  Twg  or  Twf)  until  the  process  is  completed.  The  use  of  eqs.  (8)  through 

t 

(13)  requires  a  knowledge  of  Cp,  Hch,  Hch/  MJeG/  Mch/  HV/  Hf,  Twg/  ar>d  Twf* 
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III.  EVALUATION  OF  PARAMETERS 


The  foregoing  equations  are  sufficient  to  accurately  describe  the  re-entry  behavior 
of  the  fuel  rod  up  to  the  point  where  the  rod  starts  to  break  apart.  They  are  still  valid  after 
the  break-up  begins,  provided  the  various  parameters  associated  with  the  imperfectly  shaped 
remnants  of  the  rod  can  be  calculated  or  otherwise  evaluated.  The  Assumptions  1  through  5 
are,  to  within  the  present  state  of  the  art,  quite  satisfactory.  Consequently,  the  only  items 
that  require  a  detailed  analysis  are: 

CD/  CH/  Cp' 

Hc,deG'  HdeG'  *th#  Hch,  Hf,  Hv , 

Tw#  Twg/  Twf,  c  ,  G(TW), 

MdeG;  Mj-h'  Fc ,  Fv,  and 

Assumption  6. 

It  is  not  necessarily  true  that  all  of  these  parameters  will  be  required  in  a  particular  re-entry 
computation  since  the  phenomena  are,  in  some  ways,  related  and  may  interfere  with  one 
another.  For  example,  the  formation  of  an  oxide  layer  on  the  surface  of  the  rod  may  prevent 
the  transpiration  of  the  hydrogen.  Each  of  these  parameters  will  now  be  analyzed. 

A.  VALIDITY  OF  ASSUMPTION  6,  AND  DETERMINATION  OF  Cp  AND  Tw 

It  will  now  be  shown  that  Assumption  6  is  quite  convenient  and  permits  one  to  write 
the  equations  of  heating  (eqs.  9,  11,  12,  and  13)  in  forms  which  are  readily  manageable.  By 
writing  Tw  and  Cp  in  functional  forms,  the  need  for  assuming  a  uniform  temperature  through¬ 
out  the  rod  can  be  circumvented,  as  will  be  shown  presently. 
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It  will  first  be  shown  that  the  transient  portion  of  the  heat  wave  in  the  rod  can  be 
neglected  to  within  certain  limitations.  Now,  the  major  point  of  concern  is  the  temperature 
difference  between  the  surface  and  center  of  the  rod .  Since  this  difference,  in  an  infinitely 
long  rod,  will  be  greater  than  or  equal  to  the  difference  for  a  finite  rod,  an  upper  bound  is 
obtained  by  considering  an  infinitely  long  rod.  According  to  Carslaw  and  Jaeger,8  the  tem¬ 
perature  at  a  distance  x  from  the  center  of  a  circular  cylinder  of  infinite  length,  for  ^  =  con¬ 
stant,  is  given  by 


where  TQ  is  the  initial  uniform  temperature  of  the  rod  at  the  time  the  heat  flux  qx  into  the  rod 
started  to  flow,  Jn  is  the  ordinary  Bessel  function  of  order  n  (JQ  being  the  zero-th  order),  and 
the  bn  are  the  positive  roots  of 

-Mb)  =  0.  (15) 

Here  qj  =  normal  heat  flux  into  surface,  and  pm,  Cp,  and  km,  are  respectively  the  density, 
specific  heat,  and  thermal  conductivity  of  the  rod.  The  quantity  of  interest  is  the  temperature 


difference  between  the  surface  and  center  of  the  rod,  namely. 


(16) 


The  exponential  damping  term  assures  that  the  transient  portion  of  the  temperature  difference 
is  rapidly  damped  out.  The  first  four  terms  of  the  series  are  tabulated  as  follows: 
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where 

da  =  .011ft3,  !<„,«  .007Btu/ft-sec-°R,*  380  Ibs/ft3,  Cp« 0. 1  Btu/lb -°R. 

From  this  table,  it  is  seen  that 

t[j,  *]-T(°rO-  °-13l]  for  t=  1  second, 

T[£,  t]  -  T(0,  t)  =  [  1  -  0.01 1]  for  t  =  2  seconds, 

^  =  ?  "  0.001]  for  t  =  3  seconds. 

From  this  it  is  apparent  that  the  temperature  difference  is  a  minimum  at  the  time  qx  is  "turned 
on"  and  reaches  an  approximately  equilibrium  value  in  a  few  (say  3)  seconds.  It  is  quite 
important  to  the  present  problem  that  the  transient  temperature  effect  dies  out  quickly.  For 
high  altitude  fuel  element  release  (>250,000  ft.),  no  large  changes  in  qt  over  a  period  of 
three  seconds  or  less  will  be  experienced  by  the  re-entering  rod.1*  It  is  apparent  from  the 
above  tabulation  that  the  temperature  difference  between  the  surface  and  center  of  the  rod  is 
given,  for  practical  purposes,  by 

T[i]-T(0)  -  ii»7.4q,°R  (17) 

_ *Km _ 

*  Actually,  km  varies  from  about  0. 005  Btu/ft-sec-°R  at  830°Rto  .009Btu/ft-sec-°R  at  3 1 70°R 
according  to  ARF  data. 
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wh«r«  is  the  net  heat  flux  (Btu/ft3 -sec)  into  the  sample,  that  is,  qt  =  C|_|(h«  “hw)ltoM»  + 
Hch  ”  ®  °  averaged  over  the  surface  of  the  rod . 

In  view  of  eq.  (17)  and  the  fact  that  qtis  not  expected  to  exceed  about  100  Btu/ft3  - 
sec  for  fuel  rod  release  altitudes  of  250,000  ft.  or  greater,  even  with  chemical  reactions,  it 
is  reasonable  to  assume  that  the  temperature  variation  within  the  rod  will  not  exceed  740PR. 

At  higher  temperatures  (where  1^  is  larger  than  .00 7),  this  temperature  difference  will  be  less. 
It  will  be  shown  in  the  following  treatment  of  Tw  that  either  considerably  shorter  time  incre¬ 
ments  than  three  seconds  or  greater  average  heat  fluxes  than  100  Btu/fta-sec  will  not  introduce 
excessive  errors  in  the  calculation  of  Tw. 


2.  Determination  of  Tw 

In  order  to  use  eq .  (14)  with  varying  q, ,  some  means  must  be  devised  for  accounting 
for  the  temperature  distribution  at  the  beginning  of  each  time  increment.  In  view  of  the  fact 
that  the  transient  part  of  eq .  (14)  will  die  out  in  three  seconds,  it  is  convenient  to  divide  the 
heat  flux  pulse  into  a  series  of  consecutive  heat  pulses,  q^ ,  which  appear  as  inputs  to  the  fuel 
element  at  successive  intervals  that  are  three  seconds  or  greater  apart.  Thus,  in  accordance 
with  eq.  (14),  the  wall  temperature  during  the  first  time  interval,  (At0),  is 


T  /A  *  \  -  T  4.  4<fr°Af°  4.  ilfd 

Tw<4'o)'  T°  %C.d  8k. 


(18) 


m 


Note  that  &t0  can  be  varied  continuously  from  three  seconds  onward  but  should  be  set  at  three 
seconds  up  to  the  time  At©  =  3  seconds  in  calculating  Tw  during  the  first  three  seconds.  After 
the  time  At0  has  elapsed,  a  new  heat  flux  pqlse  replaces  qlQ.  This  second  heat  pulse 
sees  a  different  value  of  the  rod  temperature  than  the  original  T0.  The  average  increase  in 
temperature  of  the  rod  during  the  first  heat  pulse  is  obtained  from  the  enthalpy  of  the  rod  at  the 
end  of  this  pulse .  If  it  is  assumed  that  the  heat  capacity  of  the  rod  remains  constant  over  the 
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temperature  range  present  in  the  rod,  the  enthalpy  of  the  rod  is 


<4 


(hr«<>total  "  2"  pmV'  i  7T(*,<f5 


Substituting  eq.  (14)  for  T(x)  into  this  expression,  holding  t  and  Cp  constant,  and  integrating 


gives, 


nd* 

OVodHotal  =  4  Pm-tCpT0  +  nd£  =  (hrod)o  +  Al\od' 


where  Ahroc|  is  the  increase  in  rod  enthalpy  due  to  the  input  qj,  for  the  time  t.  Since 

__  _  rtd3 1  — 

(^rod)^0fa|  =  Pm  CpT  *“4 —  •  'n  terms  ^e  average  rod  temperature,  T,  and  average  specific 

heat  capacity,  Cp,  it  follows  that  the  average  temperature  of  the  rod  at  the  end  of  the  first  time 

interval,  At0,  is 


T  -  T  _L  44oAto 

*1  1  o  — • 

pmCpd 

The  value  of  T$  +  l  can  be  written  immediately  as 


Ts+i  =  T°+  idL*”"1 


09) 


Here  T$  +  t  is  the  average  temperature  of  the  rod,  based  on  constant  ?p,  at  the  beginning  of 
the  time  interval  At$.  Equation  (18)  may  also  be  written  in  general  form  as 


_  »  •=  4qis  ^s  qj$d 

Tw(Ats)  ~  Ts  +  i  +  _  _  + 


RnCpd 


8  k, 


(20) 


m 


Now,  it  must  be  noted  that  Cp  varies  with  temperature.  To  take  this  into  account,  it  must  be 
considered  in  calculating  T^  +  1.  The  value  of  Cp  at  a  particular  temperature  can  be  written  as 
Bi  +  Bg T .  If  this  functional  form  is  substituted  into  eq.  (20),  the  result  is 


T/At)  =  T  +  ii»l  + -iiliili _ 

w<  *  ,+  1  ek*,  (u(M4T,+1)d 


(21) 
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The  appropriate  value  of  TJ  +  1  for  use  in  eq.  (21)  is  obtained  by  substitution  of  +  B8T,  + 
into  eq .  (19)  and  solving  for  Ts  .  The  result  is 


(22) 


When  this  value  of  TJ+1  is  substituted  into  eq.  (21),  the  wall  temperature  is  given.  This  value 
of  Tw  is  used  in  equations (8)  and  (11).  It  is  apparent  from  eq.  (22)  that  the  last  time  increment 
in  the  series  (At$)  has  a  very  small  effect  upon  Tj  +  j,,  or  on  Tw(At,),  provided  that  s  is  large 
(say  1 0  or  more) . 


3.  Determination  of  Cp<  ^  and  Ba 

In  accordance  with  the  averaging  process  represented  by  eq.  (22),  it  is  clear  that  the 
effective  value  of  Cp  is 

ep-  <»> 

This  is  valid  providing  that  At,  is  kept  sufficiently  small  and  that  the  average  value  of  Cp  con 
be  considered  constant  during  the  time  interval  At,.  Now,  it  has  been  pointed  out  that  the 
transient  effect  due  to  increasing  qx  from  q^,  to  q1$  has  a  small  effect  on  T,  for  large  values 
of  s,  that  is,  when  the  transient  heat  pulse  is  superimposed  upon  a  high  enthalpy  rod,  the  tran¬ 
sient  part  of  the  heat  pulse  only  contributes  a  small  amount  to  the  total  enthalpy.  Consequently, 
one  can  use  quite  small  values  of  At,,  certainly  on  the  order  of  one  second,  for  the  large  values 

of  s  near  the  peak  heating  portion  of  the  trajectory. 

Because  of  the  large  mass  of  the  fuel  element  relative  to  the  mass  of  the  oxide  formed 
thereon,  the  change  in  heat  capacity  due  to  the  conversion  of  some  of  the  fuel  element  to 
oxides  and/or  nitrides  can  be  neglected  to  a  good  approximation.  Hedge,  et  alx0(ARF  data) 


20 


give  the  values  of  Cp  without  the  coating  layer  of  Zr^  and/or  ZrN,  and  Wicks  and  Block30 
give  the  heat  capacities  (as  a  function  of  temperature)  for  Zr09  and  ZrN. 

It  should  be  noted  that  Hedge,  et  al,10  found  that  during  the  oxidizing  of  a  sample 
of  the  fuel  element  material  at  the  temperature  1930°R,  no  hydrogen  was  released.  It  took 
the  sample  about  160  seconds  to  come  to  an  essentially  equilibrium  temperature  in  this  particu¬ 
lar  case.  This  indicates  that  the  hydrogen  does  not  evolve  up  to  temperatures  of  1930PR.  Their 
data  on  the  oxidation  in  a  10%  air,  90%  argon  atmosphere  at  3010PR  clearly  show  an  initial 
decrease  in  mass  which  is  consistent  wirh  the  loss  of  most  of  the  hydrogen  in  the  specimen.  The 
sample  reached  equilibrium  in  about  70  seconds  in  this  case,  and  the  center  temperature  reached 
1930°R  in  about  20  seconds.  Samples  of  fuel  element  material  were  exposed  to  a  simulated  re¬ 
entry  trajectory  (by  varying  the  rate  of  heat  energy  input  simultaneously  with  atmospheric  pres¬ 
sure)  by  Littman31  at  SRI  (Stanford  Research  Institute).  He  observed  that  water  vapor  appeared 
during  these  tests.  The  GTC  chemical  and  metal lographi cal  analyses  of  the  SRI  samples  confirm 
the  loss  of  hydrogen  from  the  samples  that  underwent  the  simulated  re-entry  tests.  Except  for 
the  very  low  heating  trajectory  (SRI  Trajectory  1A)  simulated  by  SRI,  essentially  all  of  the  hy¬ 
drogen  escaped,  and  even  on  this  low  heating  trajectory,  about  75%  of  the  hydrogen  escaped. 
These  data  represent  the  only  experimental  data  available  with  which  to  determine  the  rate  at 
which  the  hydrogen  escapes  through  the  surface.  Although  ARF  tried  to  establish  the  diffusion 
rate  of  hydrogen  through  an  oxide  coating  on  a  sample  specimen  of  fuel  element  material,  the 
experiments  failed  because  of  loss  of  oxide  through  structural  failures.  This  wiH  be  discussed 
in  more  detail  in  the  next  subsection  under  G(TW) . 

There  are  two  ways  of  handling  the  specific  heat  capacity  of  the  fuel  element  material. 
The  oxides  and/or  nitrides,  if  shown  to  be  significant,  can  be  handled  separately.  The  ARF 
data  show  that  the  heat  capacity  of  the  hydrided  material  is  higher  than  that  of  the  dehydrided 
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material  at  temperatures  above  about  700^R  and  equal  to  it  at  high  temperatures  (2040PR) . 

It  is  here  assumed  that  no  hydrogen  is  lost  during  re-entry  at  temperatures  below  2040PR. 

The  above  mentioned  data,  as  well  as  the  analytical  evaluation  of  such  data  (see  Appendix  II) 

support  this  position .  Thus,  at  high  temperatures,  it  can  be  shown  from  the  ARF  data  that  the 

heat  capacity  of  the  dehydrided  material  is 

Cp|  =  0.0433  4  1 .33  x  10_ST  Itu/lb  -  °R  for  T  ^  2040°R.  (24) 

!Zr:U 

Hence,  Bi  =  .0433,  and  Bg  =  1 .33  x  10“sin  this  temperature  range.  This  will  be  slightly  in 
error  in  the  vicinity  of  the  transition  point  (a  to  S  phase  at  2040°R)  but  will  be  accurate  at 
the  high  temperatures.  Similarly,  the  heat  capacity  of  the  hydrided  material  is 

=  0.108  4  2.65  x  10"5T  Btu/lb  -  °R  for  700  <T<2040°R .  (25) 

Zr:H:U 

Hence,  Bx  =  0.108,  and  Bg  =  2.65  x  10”B  in  this  temperature  range. 

Since  the  value  of  Cp  does  not  change  markedly  in  going  from  eq.  (24)  to  (25),  it  is  a 
fair  approximation  to  use  either  of  these  values  in  eqs.  (21)  and  (22).  For  temperature  very  near 
to  2040°R/  eq.  (14)  can  be  used  to  establish  what  portion  of  the  fuel  rod  should  be  treated  as 
having  a  heat  capacity  given  by  eq.  (24),  with  the  remainder  being  treated  as  having  a  heat 
capacity  given  by  eq .  (25) . 


B.  DRAG  COEFFICIENT  AND  VARIATION  WITH  ANGLE  OF  ATTACK 

The  drag  of  the  fuel  element  rod  can  be  calculated  from  the  "modified"  Newtonian 
theory.  For  this  calculation,  it  is  convenient  to  break  the  total  drag  down  into  two  parts;  that 
is,  the  part  due  to  the  radial  cylindrical  surf6ce  and  the  part  due  to  the  ends. 

Raymond's  drag  curve8  is  of  little  value  in  the  present  application  since  the  i/6  ratio 
is  quite  different  from  that  for  the  fuel  element  rod.  Fortunately,  Penland11  has  the  Newtonian 
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Cq  curve  for  an  infinite  cylinder.  This  curve  i»  plotted  in  Figure  3.  For  a  finite  cylinder  with 
flat  or  hemispherical  ends,  the  values  must  be  modified  to  allow  for  the  end  effects.  The  local 
pressure  on  the  end  face  of  the  rod,  which  !s  now  treated  as  a  flat  disc,  is  equal  to  Pq  cos *0, 
where  P^  is  the  normal  shock  recovery  pressure .  Thus,  the  normal  force  =  l&  cos*0A<j,  where 
Aj  is  the  area  of  the  disc .  The  C|)j  for  the  disc,  based  on  the  area  of  the  disc,  is  given  by 


P'  cos COS0  Aj  Pq  cos*0 

cDd  -  j — - r - 7  • 

S  do  3 

Here  is  the  free  stream  dynamic  pressure.  The  drag  coefficient,  based  on  the  broad¬ 

side  area  A  of  the  cylinder  is  given  by  the  relation 

cDd  =  , 


cDd 


CDd 


TT  1  .29  P, 

4  x 


Y  p 
2P® 


oft  l 


<3 


/ 


CDd=  .Uicm’t, 

wh*r*  Z  *  TS  ■ 

Table  1  shows  how  the  total  drag  coefficient  depends  on  the  pitch  angle,  0.  The  drag 
coefficient  for  the  end  disc  is  also  plotted  in  Figure  3.  Finally,  the  total  drag  coefficient  is 
shown  in  Figure  3.  A  good  analytical  fit  to  this  total  drag  curve  is  given  by  the  following 
expression: 

CD  tota|  =  .144  cos30  +  1 .2  slna0  -  .175  sin*(20).  (26) 

This  function  is  tabulated  in  Table  2  and  plotted  in  Figure  3. 
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Based  or  modified  Newtonian  flow,  see  Ref.  11  . 


o 


Angle  of  Attack  0 
Figure  3.  Cylinder  Drag  Coefficient 


Table  1 .  Components  and  Total  Drag  Coefficients  of  a  Cylinder  with  1/6  *  10 


-a* 


0 

COS  0 

COS®0 

Cq  Disc 

Cq  Cylinder 

Cq  Total 

QO 

1 

1 

.146 

0 

.146 

10 

.9848 

.94 

.137 

.01 

.14? 

20 

.9397 

.82 

.120 

.05 

.170 

30 

.8660 

.643 

.0938 

.155 

.249 

40 

.7660 

.446 

.0652 

.34 

.41 

50 

.6428 

.263 

.0384 

.54 

.58 

60 

.5000 

.125 

.0183 

.78 

.80 

70 

.3420 

.0400 

.00584 

1.0 

1.0 

.1736 

.0052 

.00076 

1.15 

1.15 

90 

0 

0 

0 

1.2 

1.2 

Table  2.  Tabulation  of  Equation  26 


<t> 

.146  cos*0 

1 .2  s»na  0 

.175sin^2fl) 

Cq  Total 

0 

.146 

0 

0 

.146 

10 

.137 

.036 

.0205 

.1525 

20 

.120 

.14 

.0723 

.1877 

30 

.0938 

.30 

.131 

.2628 

.0652 

.496 

.169 

.3922 

50 

.0384 

.704 

.169 

.5734 

60 

.0183 

.90 

.131 

.7873 

70 

.00584 

1.058 

.0723 

.992 

80 

.00076 

1.16 

.0205 

1.14 

90 

0 

1.20 

0 

1.20 

If  the  rod  has  hemispherical  ends,  then  the  ends  have  a  constant  drag  coefficient.  At 
hypersonic  speeds,  this  value,  according  to  Kuehn,4*  Is  Cq  ■  0.94  if  based  on  cross  sectional 
area  of  the  sphere .  It  should  be  noted  that  if  the  cylinder  broadside  area  Is  used,  then  a  dif¬ 
ferent  area  from  that  used  in  the  last  analysis  should  be  selected.  The  new  area  would  simply 

be  (l  -  d)d  +  ^r-  .  The  drag  coefficient  for  this  case  Is  obtained  by  multiplying  numbers  in 

the  fifth  column  of  Table  1  by  the  ratio  £(6  -  d)d  +  l^>d .  This  modifies  eq .  (26)  toread 
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(27) 


CD  total  I  =  0.069  4  1.17  sin8  0  -  0.171  si na (2 0) . 

■spherical 
caps 

Thus,  if  or  when  the  edges  of  the  fuel  element  burn  away,  eq.  (27)  is  more  appropriate  than 
eq .  (26) . 

C.  EVALUATION  OF  HEAT  TRANSFER  COEFFICIENT  Ch 

In  the  re-entry  regime  of  interest  for  the  fuel  element  rod  (altitudes  of  300,000  feet  or 
less),  the  heat  transfer  can  be  considered  to  be  due  to  laminar  aerodynamic  heating.  The  heat 
transfer  coefficient  for  the  "no-blowing"  case  can  be  derived  from  the  report  of  Vidal  and 
Wittliff.4  The  effects  of  blowing  are  then  handled  in  accordance  with  eqs.  (9)  through  (13). 
The  empirical  expressions  derived  closely  approximate  the  experimental  results  obtained  in  the 
Cornell  Aeronautical  Laboratory^  Six  Foot  Hypersonic  Shock  Tunnel .  These  tests  were  carried 
out  at  Mach  number  from  17.13  to  23.46  at  stagnation  temperatures  varying  from  1800®K  to 
5100^K  for  a  wide  range  of  Reynolds  numbers.  These  Reynolds  numbers  cover  those  appropriate 
to  the  fuel  element  over  the  altitude  range  150,000  feet  to  300,000  feet.  The  experimental 
data  include  both  the  real  gas  effects  and  the  hypersonic  interaction  effects  which  will  be  en¬ 
countered  by  the  re-entering  fuel  element.  The  chemical  effects  of  reaction  heating,  trans¬ 
piration  cooling,  etc.,  will  have  to  be  handled  independently.  The  determination  of  the 
aerodynamic  heating  effects  is  established  in  the  following  six  subparagraphs. 


1 .  Definition  of  Heat  Transfer  Coefficient 

The  laminar  heat  transfer  coefficient  is  defined  as19 


CH(K8,*,  *) 


. i.  : 

<boV®  (h*  "  K#) 


(28) 
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This  empirical  expression  was  derived  by  fitting  a  cubic  equation  to  the  data.  Equation  (35) 
is  quite  valid  in  the  range  0.1  ^  K8  ^  100.  It  agrees  with  the  viscous  shock  layer  theory  of 
Cheng  as  well  as  that  of  Cohen  and  Reshotko.4 

It  should  be  noted  that  for  K^O.I ,  one  gets  into  the  transition  and  free  molecular 
flow  regimes/  and  eq.  (35)  is  not  applicable.  Also,  for  very  large  values  of  K8,  eq.  (35)  fails. 
It  is  actually  theoretically  correct  for  the  conditions:  Tj3>Tw,  Pr  =  0.71 ,  and  e*  =  0.10; 
where  &  is  the  density  ratio  across  the  normal  shock,  and  Pr  is  the  Prandtl  number.  However, 
it  does  fit  the  experimental  data  quite  well,  as  indicated  by  Figure  4. 

3.  Heat  Transfer  Distribution  Around  the  Circumference  of  a  Cylinder 

In  the  following  discussion,  the  position  around  the  circumference  of  the  rod  will  be 
designated  by  6,  with  6  -  0  at  the  stagnation  point  for  the  rod  in  cross  flow.  The  quantity 
q(0)  will  be  used  to  denote  the  heat  transfer  rate  to  the  cylindrical  surface  at  the  angle  8  . 

It  should  be  pointed  out  that  there  are  discrepancies  between  the  various  experimental 
values  of  the  ratio  of  q(§")/q$  for  various  large  values  of  ?  and  some  discrepancies  for  various 
values  of  the  pitch  angles  0 .  However,  it  can  be  shown  that  the  scatter  in  the  experimental 
data  on  variations  with  $  are  submerged  in  the  effects  of  variations  in  the  data  on  K8 .  Con¬ 
sequently,  one  is  justified  in  ignoring  the  effects  of  0  on  q(8)  until  better  experimental  data 
become  available.  The  arguments  for  using  a  single  parameter  function  for  q(  8)/qs  or*  a* 
follows: 

a .  The  laminar  heat  transfer  distribution  around  the  circumference  of  a  cylinder  re¬ 
entering  with  its  cylindrical  axis  normal  to  the  flow  (0“  90°)  is  shown  in  Figure  5  (see  also 
Figure  6,  Ref.  4),  and  results  for  0  up  to  45.5°  are  shown  in  Reference  4,  Figure  8.  The 
experimental  values  obtained  by  Vidal  and  Wittllff  (for  M*  17  -  23)  ore  considerably  lower 
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Figure  4.  Heat,  transfer  to  a  yawed  cylinder 


than  those  obtained  by  Tewflk  and  Giedt.16  This  Is  particularly  true  for  the  back  side  of  the 
rod,  i.e.,  for  90°<?<180°. 

b.  For  0<  0  <90°,  the  Vidal  and  Wittliff  data  agree  well  with  the  theoretical  distribu¬ 
tion  derived  by  Beckwith10  for  laminar  boundary  heating  of  a  cylinder.  (See  Ref.  4  for  detailed 
comparison.) 

c .  Laminar  separation  on  the  back  side  of  the  cylinder  is  evidenced  by  the  dependence 
of  q(9)on  K*  (see  Fig.  5).  That  is,  laminar  separation  is  affected  by  the  Reynolds  number. 
Fortunately,  the  effect  of  laminar  separation  is  only  pronounced  for  0  >  120°  for  which 
q(9)/q,<.03,  l.e.,  less  than  3%  of  the  stagnation  value.  Thus,  one  is  justified  in  fitting  an 
average  curve  through  all  of  the  experimental  data  points  for  various  values  of  Ka  . 

d  .  Although  the  theory  of  Beckwith  predicts  a  dependence  of  q(  @)  on  0  ,  the  data  of 
Vidal  and  Wittliff  (see  Figure  5)  shows  that  if  there  is  such  a  dependence,  it  is  submerged  by 
the  influence  of  laminar  separation.  The  spread  in  data  is  represented  by  the  fact  that  for 
0  =  90°  and  45 .5°,  observed  values  of  q(0Vq*  are  essentially  the  same  for  values  of  K*=  2 
and  5,  whereas  at  0  =  60°  and  K8  =  1.6,  the  observed  value  is  quite  different.  Hence,  it 
may  be  concluded  that  variations  in  q(0)/q$  with  0  are  not  justifiable  on  the  basis  of  experi¬ 
mental  data . 

e.  Finally,  in  view  of  arguments  (c)  and  (d),  it  is  concluded  that  a  single  parameter 
curve  of  q($)/qt  versus  0  is  all  that  can  be  justifiably  employed .  A  curve  that  fits  the  data 
presented  by  Vidal  and  Wittliff4  is 

•°9io[q(^/9s]  *  3[cos|-l]  .  (36) 

This  relation  is  compared  with  the  aforementioned  experimental  results  in  Fig.  5.  It  is  seen 
that  the  empirical  equation  agrees  quite  well  with  the  experimental  data  for  the  forward 
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portion  of  the  fuel  rod  and  is  in  fair  agreement  on  the  back  side  where  the  heat  transfer  is 
quite  small  compared  to  qs. 


It  will  be  assumed,  as  mentioned  above,  that  the  fuel  rod  will  spin  about  its  cylindri' 
cal  axis,  hence  it  is  practical  to  use  an  average  value,  q„,  for  the  heat  flux.  Rewriting 
eq.  (36)  in  exponential  form  and  integrating  gives 

qa  =  £*  ^[6.9078  (cos  a -1)]  da  =  0.310  qs.  (37) 

This  integral  was  evaluated  from  Figure  6,  which  illustrates  the  variation  of  q  with  0\ 


4.  Heat  Transfer  Rate  to  a  Pitching  Cylinder 

Now  that  the  heat  transfer  at  the  stagnation  point  (eq.  35)  and  the  effect  of  spin 
(eq .  37)  on  it  have  been  developed  in  analytical  form,  there  remains  to  be  determined  the 
effect  of  pitch  angle  on  the  heat  transfer. 

The  experimental  data  on  stagnation  line  heat  transfer  to  yawed  cylinders  have  been 
reviewed  in  Ref.  4  and  presented  in  Figure  7  thereof.  This  data,  together  with  the  Vidal- 
Wittliff  data  was  used  to  construct  the  curve  shown  in  Figure  4.  Theoretical  solutions  have 
been  obtained  by  Cheng  and  Chang17  which  indicate  (C^t)/sin0  is  principally  a  function  of 
K3^),  where 

R#oo 

**(6)  =  -  .  (38) 

Y00M800C*U)|»in*| 

Similarity  with  equations  (31),  (32),  and  (33)  gives 


r*(^\  =  m*(*)Tco 

(39) 

umT*  (») 

T *(«)  -  i(Tw  +  T,.ln*«), 

(4Q) 
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and  n*(0)  is  to  be  evaluated  at  the  reference  temperature ,  T*(0).  The  sin8  0  term  arises  from 
viscous  shock  wave  theory.17  In  terms  of  more  conventional  parameters,  eq.  (38)  may  be 


written  as 


[_Tw  +  T$sin80] 

- 1 

jc 

*5! 

J 

for  0.1  ^  K8(0)  ^  100. 


(41) 


The  absolute  value  is  used  on  the  sin0  in  order  that  0  may  vary  from  0  to  2n  .  The  limits 
(0.1  to  100)  for  the  range  of  Ks(0)  are  in  good  agreement  with  the  experimental  data  of  Vidal 
and  Wittliff .  The  curve  must  be  extrapolated  beyond  the  upper  limit  to  K8(0)  =  100. 

Although  Cheng  and  Chang  found  that  C^y/sin0  is  weakly  dependent  upon  0,  both 
theory  and  the  experimental  results  of  Vidal  and  Wittliff  indicate  that  it  is  small  enough  to  be 
neglected.  This  is  referred  to  as  the  "yaw  -  independence  principle."  In  view  of  this  principle, 
eq.  (35)  still  holds  when  K®(0)  and  C^|t/sin0  are  substituted  for  C^$  and  K8,  that  is 


l°9ib(CHs/s‘n0)  =  “•337 -  .4063  logw K#(0)-  .  1 235 (jog10 K8(0)J  +  .04283  ^0 K8(0)]3 

for  0.1  ^  K8(0)^  100 K8(0).  (42) 


For  values  of  K*(0)  >100,  the  following  equation  should  be  used. 


loglo(CHj/*'n0)  =  -0.5  log  K8(0)- 0.30.  (43) 

This  equation  is  illustrated  by  the  section  of  the  curve  in  Figure  4  that  lies  between  K8^)  = 
100 and  1000.  For  K®(0)<O.l,  the  value  of  C^*  calculated  at  K8(0)=  0.1  by  eq.  (42) 
should  be  used,  i  .e . , 

CHs(K3/*)  =  CH,(O.1,0).  (44) 

Equation  (44)  is  not  needed  for  altitudes  below  300,000  ft. 
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5.  Heot  Transfer  to  Afterbody 

For  values  of  K*(0)  greater  than  some  value  l?*(0),  C^,  Is  equal  to  the  afterbody 
heat  transfer  coefficient  (neglecting  ends  of  cylinder)  of  the  cylinder  in  the  end-on  flow 
position.  (The  heat  transfer  to  the  ends  of  the  cylinder  will  be  treated  in  the  following  sec¬ 
tion.)  This  is  a  little  more  involved  and  requires  some  detailed  discussion.  Based  on  Blasius 
incompressible  flow  skin  friction,  modified  for  application  to  compressible  flow  by  evaluating 
at  the  reference  enthalpy  of  Eckert, 18  together  with  Colburn's  version  of  the  Reynold's 
analogy,19 

S,  *  -j-  Cf  Pr'*’, 

qab“  ch».vA-M-  MS 

Here  q^  =  heat  flux  to  afterbody, 

|i  =  viscosity  of  air, 

Pr  =  Prandt I  number  *  4^9  Y  -5), 

Cf  =  local  skin  friction  coefficient  (Cf  =  mean  value  of  Cf),  and 
hr  -  recovery  enthalpy  (defined  below) . 

The  value  of  CH  •*  given  by  Hayes  and  Probstein  as19 


0.664  -»/,  /(PH)*" 

:Hab=^Pr  ' 


<*) 


The  asterisk  refers  to  quantities  evaluated  at  the  reference  enthalpy,  the  subscript  e  refers  to 
local  values  at  the  outer  edge  of  the  boundary  layer,  and 


<*•>*  = 


w 

^e 


(47) 


The  reference  enthalpy  of  Eckert  is  used,  i  .e. , 
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h*  =  ^e4  h*)  + 0.22(1* -h*). 


(48) 


and 

hr  =  he  +  <hs  -  he)  >/pT  •  (49) 


Cohen80  gives  the  value  of  ( p|i)V(PM)e  as 

(PM)*  1  -  1 .0213 

1  -  (he/h1)0**'89] 

Me 

1  -  1.0213 

[l  -  (hVh1)0*3889] 

The  value  of  h'  is  272,234  Btu/slug  (8465  Btu/lb). 

Now  the  values  of  ,  Ve,  Me,  and  h^  (local  flow  on  the  afterbody)  must  be  found. 
This  can  be  accomplished  as  follows:  For  present  purposes,  perfect  gas  theory  will  be  used,81 
and  the  pressure,  Ps,  behind  the  shock  wave  (edge  of  boundary  layer  in  front  of  stagnation 
point)  will  be  calculated  from  the  Rayleigh-Pitot  tube  formula,  which  gives 


(51) 


The  value  of  ,  the  local  Pitot  pressure  at  the  outer  edge  of  the  boundary  layer,  can  be  deter¬ 
mined  from  the  graph  presented  by  Chernyi,88  which  is  based  on  hypersonic  flow.  An  analytical 
fit  to  this  graph  is 


k'(y) 


(52) 


Here  x0  is  the  distance  along  the  surface  of  the  cylinder  from  the  stagnation  point.  The  quantity 
K'(y)  =  0.095  for  y  =1-4;  =  0.080  for  y  *  1 .2;  and  =  0.066  for  Y“  1  -1  •  A  reasonable  ap¬ 
proximation  is  to  take  d/xe  =  .2  (the  average  value),  y*  1*2,  and  CD  3  1  *84  (see  ref.  43), 


I 


Pa  =  + 2.02x10-^^. 


(53a) 


It  should  b«  mentioned  that  Kuehn  «  used  Pa/2  in  place  of  P®  in  this  last  equation .  For 
M00>  15,  this  makes  little  difference.  For  smaller  values  of  M®,  one  should  use 

Pe  =  rPa  +  2.02x  lO-3^^  .  {i 

For  an  isentropic  expansion  of  the  Pitot  pressure  along  the  afterbody. 


The  temperature,  Te,  in  terms  of  the  temperature,  Ts  (at  the  edge  of  the  boundary  layer  in  front 
of  the  stagnation  pcint),  Is  given  by 


t.  r,  +  i±T*  rrei3/7 
js  L  5  J  /  °r  Te  -  Ts  [_  J  , 

where  Ts  *  h/Cpg,  and  Cpa  is  the  specific  heat 
sound  speed  is  ae,  and  Ve  =  <>eMe/  where 


fe  -\aA 


capacity  of  air  =7.73  Btu/slug  -  °R . 


(55) 

The 


ae  =VYRgV  •  Rg  =  1715  ft-lbs/slug  -°R.  (53) 

Therefore,  from  the  above  relations,  it  is  readily  shown  that 

Ve  =  V^YRg(h$Cpa  “  1)  •  (57) 

The  density  is  Pe  =  P#/RgT#,  which,  by  the  above  formulas,  is 

0  .  6M«P®  rif®  +  0J26pa>V2o  1BA 

#  SRgh.Cpo  L  7M»oo-l  J  •  (58) 

The  viscosity,  according  to  Sutherland's  formula,  becomes 
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Here,  the  values  of  C  and  S*  have  been  taken  from  Reference  21  (see  eq.  33).  The  value  of 


h«  Is  taken  as  he  =  Te^pa/  that  Is 


When  equations  (57)  through  (60)  are  substituted  into  (46)  through  (50),  all  of  the  parameter* 
in  eq.  (45)  are  given  in  terms  involving  conditions  in  front  of  the  shock  wave  and  at  the  stag¬ 


nation  point  enthalpy.  This  latter  quantity  is  simply 


L  _  ^0  —  C  T  (61) 

hs  “  U  CPaT*  * 

Note,  now,  that  C^b  must  be  used  in  place  of  the  value  of  Chs/  a*  calculated  from  eq.  (33), 
when  the  latter  becomes  small .  That  is,  use  C^b  when 

CHab-  O.SIOCh.  •  m 


6.  Heat  Transfer  to  Ends  of  Cylinder 

The  heat  transfer  rate  to  the  end  face  of  a  cylinder  in  end-on  flow  is  nearly  constant 
according  to  Kemp,  Rose,  and  Detra  (see  Fig.  5  of  Ref.  23).  According  to  Detra  and  Hidalgo/ 


the  heating  rate  is 
865 

qend  *  (d/2) 


V,  k^n  [  rr]  flow)  • 


Here  C  *  (du/dx),  *  stagnation  point  velocity  gradient, 


£n  «  Newtonian  value  of  C  /  ond 

*  wall  gas  enthalpy  at  Tw  ■  300®K  (540®R),  hence,  hyy^  ■  CpgT  m  7.73x  540  * 
W  4175  Btu/slug. 
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This  equation  is  accurate  to  within  10%  for  6000^  26,000  ft/sec  and  for  8  x  10“*  <6 


For  large  values  of  M,  Probstein86  shows  that  the  correct  value  of  C  i* 


1.2 


C  =  "3“ 


Too 


(65) 


Hence , 


(66) 


Since  the  heat  transfer  rate  is  expected  to  be  proportional  to  Ps,  and  the  pressure  distribution 
on  the  end  of  the  cylinder  varies  as  cos8  0  (see  Weinstein48  for  experimental  verification), 
the  heating  rate  becomes,  upon  substitution  of  eqs.  (64),  (65),  and  (66)  into  (63), 


®  Hr*  **  ® 


(67) 


The  appropriate  values  of  Y  have  been  tabulated  by  Wittliff  and  Curtis;86  however,  Y  *  1 .2 
would  not  lead  to  gross  error  in  the  high  heating  part  of  the  re-entry  trajectory. 

Summarizing,  the  values  of  CH  to  be  used  in  equations  (8),  (9),  and  (11)  are: 


a.  0.1  ^  K3(0)£  100<K\ 

Take  the  sum  of  the  value  obtained  by  using  eq. (42)  (multiplied  by  0.310) 
and  the  value  qend(^/pbo^oo^  Aend(^)  **  calculated  from 

eq.  (67). 


b.  K?(0)<O.l. 

Take  the  value  obtained  when  K8(0)  “0.1  is  substituted  into  eq.  (42),  multi¬ 
ply  by  0.310,  add  the  value  of  qend(^)/lbo^co^*  ”  ^  obtained  from  eq.(67). 
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C.  K®(0)>  100  and  Cho^  0.310  CH|. 

Multiply  the  value  of  CHs  obtained  from  eq.  (43)  by  0.310,  and  add  it  to 
W^/PcoVhj-M. 
d.  CHab  ^  0.310 CHs. 

Take  the  value  of  C^ab  obtained  from  eq.  (46)  and  add  it  to 
4end^)/ Rgo  ^ao  (^$  ~  ^w)  • 

D.  EVALUATiON  OF  e 

The  data  on  the  emittance  from  Hedge,  et  al  can  be  represented  to  within  about  1% 
by  the  equation 

e  =  -0.08+  / 2  x  10*  +  3500T  -  T8  .  (68) 

E.  EVALUATION  OF  Hch,  Hcb,  AND^|  ^ 

In  this  category,  there  are  three  major  considerations,  that  of  the  heat  input  due  to 
the  oxidation  and/or  nitridation  of  the  fuel  rod  material,  the  rate  of  weight  accretion  due  to 
these  reactions,  and  the  heat  input  (or  loss)  due  to  the  dehydriding  and  combustion  of  hydrogen 
contained  in  the  fuel  element  material.  The  available  data  have  been  analyzed  in  Appendix  II. 

,  c  .  .  ,  dWI 

I.  Evaluation  of  -rr  | 

-  dt  |ch 

The  data10 on  the  oxidation  and/or  nitridation  of  the  fuel  element  material  shows  that 
the  weight  gain  rate  during  oxidation  follows  the  parabolic  law,  except  during  dehydriding, 
within  the  accuracy  of  the  experiments.  (See  Appendix  II  for  a  detailed  analysis  of  the  avail- 

(69) 


able  data.)  That  is, 
(6m,)a=  kjt 
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where  Am,  i*  the  weight  gain  per  unit  area  in  a  time  t,  and  k,  it  a  temperature  dependent 
parameter  (reaction  constant)  which  obeys  the  Arrhenius  relation 

k,-A,.*Q,/T.  (70) 

The  oxidation  data  of  reference  10,  after  correcting  for  sample  shape  and  size,  are  plotted 
in  Figure  7.  The  straight  line  represents  the  parabolic  curve,  and  the  points  represent  the 
experimental  data.  The  values  of  A,  and  Qs  were  determined  by  plotting  logk  versus  1/T  as 
shown  in  Figure  8.  This  figure  shows  that  the  pre-factor,  A$,and  the  activation  energy,  Q$ 
(actually  RgQs),  change  at  the  phase  change  of  the  fuel  material,  i  .e.,  at  about  2040°R. 
For  t  in  minutes  and  Ams  in  grams  per  square  centimeter,  the  numerical  values  of  A,  and  Qt 


As  =  0.356  gmVcm4 -min,  and  Q,  =  2.59  x  104OR  for  T^  2040^R, 
As  *  92.5  grr^/cm4 -  min,  and  Qs  =  3.79  x  104OR  for  T  ^  2040°R. 


(71a) 


(71b) 


To  convert  A  m,  to  pounds  per  sq.  ft.,  one  first  divides  A$  by  sixty  to  obtain  grams8/ 
cm4 -sec  and  then  multiplies  by  2.04  to  convert  to  lbs*/ft4—  sec.  Thus,  eq.  (69)  becomes 

Ann  *  3.4  x  10’*(k|t)3^  Ibs/ft8,  for  t  in  seconds.  (72) 

Differentiation  of  eq.  (72)  gives  the  rate  of  weight  gain  due  to  the  parabolic  oxidation  and/or 
nitridati on  accretion,  as 


■=  -  1.7x  10-*(k/t)^.  (73) 

This  expression  is  to  be  multiplied  by  the  surface  area  and  then  substituted  into  the  re-entry 
equations  for  the  accretion  component  of  dW/dt,  i.e.,  equations  (2),  (7),  and  (12a). 
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Figure  7,  Effect  of  temperature  on  the  oxidation  of  zirconium-uranium  hydride  in  100%  air. 
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Y  x  10  4  (T  in  degrees  Rankine) 

Figure  8.  Variation  with  temperature  of  reaction  rate  constant  for  oxidotion  of  zirconium- 
uranium  hydride. 
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2.  Evaluation  of  Hch 


For  the  reaction  Zr  +  20  — ZrO^ ,  the  heat  energy  released  is,  according  to 
Wicks  and  Block,90  for  T  in  degrees  Kelvin,  given  by  the  following  equations. 

HZrO*  =  262,960  -  2.65T  +  1 .6  x  10~* T®  -  2.09  x  1(f  T^ca l/mole  of  ZrO* 

for  298°K  —  T  —  1 , 135°K .  (74) 

HZrO^  =  264,360  -  2.21  T  -  4  x  10”4Ta  -2.96  x  lOPT"1 ,  1,135°K4T  4  1,478°K.  (75) 

HZrQ8=  262,400  -3.37T  + 5.0  xlO-4T3+4. Ox  104rl,  l ,478°K4T 4  2, 100°K.  (76) 

Wicks  and  Block  also  give  the  following  relations  for  the  heat  released  in  the  reaction 
Zr  +  N  — ►  ZrN . 

HZrN  =  87,870  *  °.94T  -  2.5xlO-*T3-8.5x  1°4T “x  ,  298°K  4  T  <  1,135°K.  (77) 

HzrN  "  89, 1 10  -  0.40T  -  5.8  x  lO^T8- 1 .72x  KfT"'’1 ,  1 ,135°K  4  T  4  170GPK.  (78) 

These  values  are  in  calories  per  mole,  with  T  in  °K.  The  equivalent  increase  in  weight  due  to 
the  formation  of  a  mole  of  ZrO^  over  that  due  to  the  formation  of  a  mole  of  ZrN  is  32/14  = 
2.28.  When  this  last  expression  is  multiplied  by  2.28,  the  result  is 

2.28  HZrN  =  203,000  -  0.915  T  -  1 .33  x  10“*Ta-  3.93  x  lOV1  . 

These  relations  show  that  the  heat  energy  input  to  the  fuel  rod  per  unit  weight  gain  is  not 
markedly  different  for  the  formation  of  ZrO^  and  ZrN.  Further,  the  heat  of  formation  of 
UQa  is  almost  identical  with  that  of  ZrOB .  Therefore,  it  makes  little  difference  which  re¬ 
action  is  responsible  for  the  weight  gain  insofar  as  Hch  is  concerned .  The  ratio  of  the  heat 
energies  released  per  unit  weight  gain  by  the  fuel  element,  thatis,(HZro/2.28  HZr|vj)  for 
T  =  300,  600,  900,  1200,  1500,  and  1700°K, are,  respectively,  1.32,  1.32,  1.32,  1.31, 
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1 .31 ,  and  1 .31 .  Further/  chemical  analysis  on  samples  of  fuel  element  material  (see 
Appendices  I  and  II  show  that/  at  most,  a  few  percent  of  the  weight  gain  can  be  attributed 
to  ZrN.  Thus,  the  use  of  only  the  heats  of  reaction  of  ZrOa  will  result  in  errors  of  less  than 
a  percent  or  so.  Thar  is,  if  all  of  the  reactions  are  considered  to  lead  to  the  formation  of 
ZrO*  /  the  error  would  be  less  than  about  1%. 


3.  Determination  of  Hc^ 

The  value  of  Hcj,  is  determined  by  eq.  (12a)/  together  with  the  values  of  H^,  M,., 
and  That  is  / 


H 


H 


ch 


£h_ 


dWt 


g(H-Mch)  * 


(79) 


ch 


The  value  of  M,.  is  simply  the  molecular  mass  in  slugs  of  Zr,  i.e./  91 .22/454  x  32.2=  6x  10"? 
slugs/mole,  while  the  appropriate  value  of  is  ^ZrOg  •  0/1  123.33/454  x  32.2  =  8.45  x  lO-' 
slugs/ mole .  The  value  of  si  is  given  by  eq.  (73),  and  Hch  is  given  by  eqs.  (74),  (75),  and 


(76). 

F.  DETERMINATION  OF  G(TW) 


This  parameter,  which  constitutes  the  rate  of  loss  of  hydrogen  from  the  fuel  element, 
represents  one  of  the  weakest  theoretical  links  in  the  re-entry  model  in  that  the  data  available 
were  insufficient  to  evaluate  the  parameters  in  G(TW)  with  any  high  degree  of  confidence. 

An  insight  into  the  magnitude  of  the  effects  due  to  transpiration  cooling  can  be  gained  as  fol¬ 
lows.  If  all  of  tH*  hydrogen  is  presumed  to  transpire  during  the  peak  heating  portion  of  the  re¬ 
entry  tra(ectory,  then  the  calculated  net  heat  input  to  the  rod  would  be  less  than  that  in  the 
actual  re-entry  case,  whereas,  if  one  presumes  that  no  hydrogen  is  released  during  the  re-entry. 
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the  calculated  net  heat  input  to  the  rod  would  be  greater  than  that  in  an  actual  re-entry  case. 
An  elementary  calculation,  based  on  eq.  (9)  and  the  dynamic  heating  curves  for  a  fuel  element 
presented  in  the  North  American  Aviation  Corporation's  report  (ref.  33)  shows  that  these  two 
bounjs  differ  markedly.  That  is,  if  the  entire  hydrogen  content  of  the  fuel  element  were  to 
transpire  during  the  peak  heating  portion  of  the  re-entry  trajectory,  the  maximum  temperature 
of  the  fuel  element  would  be  about  3000°R,  whereas,  without  hydrogen  transpiration,  the  fuel 
element  is  expected  to  reach  its  melting  temperature  (calculations  based  on  re-entry  from 
300,000  feet,  1G00°R  initial  fuel  rod  temperature,  re-entry  angle  of  -0.25  degrees,  and  a 
re-entry  velocity  of  24,200  ft/sec) .  Now  that  the  possible  magnitude  of  this  effect  has  been 
put  in  proper  perspective,  the  mathematical  representation  for  the  rate  of  hydrogen  loss  will 
be  presented . 

The  process  of  hydrogen  loss  is  considered  to  be  the  following  one:  The  hydrogen  is 
released  from  the  fuel1  material  in  accordance  with  a  linear  rate  constant  X  which  depends 
upon  the  temperature  of  the  material .  The  amount  of  hydrogen  released  per  unit  area  of  fuel 
element  surface  in  a  time  t  is 

MH,  ’  Mj,  -«~U]  (80) 

where  M^0  is  the  initial  hydrogen  mass  per  unit  area  of  fuel  element  material .  (See  Appendix  II 
for  a  complete  formulation  of  this  problem.)  The  ARF  data  provides  a  means  whereby  X  can  be 
evaluated.  These  data  shew  that 

X  «  1 .01  x  I0"*T  -  1 .92,  for  T  >  ]900°R,  and  t  in  minutes.  (81) 

The  value  of  X  is  taken  as  zero  for  T <1900°R.  In  the  meantime,  an  oxide  (and/or  nitride) 
layer  has  been  building  up  on  the  surface  of  the  rod,  and  the  hydrogen  must  diffuse  through  this 
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layer  in  order  to  escape.  This  gives  rise  to  a  build-up  of  released  hydrogen  concentratin' , 

M^r  -  Am|^(t),  within  the  boundaries  of  the  oxide  layer.  Here,  A m^(t)  is  the  mass  of  the 
hydrogen  (per  unit  area  of  surface)  that  has  escaped  through  the  surface.  Since  the  concen¬ 
tration  of  hydrogen  outside  of  the  oxide  layer  can  be  considered  as  being  essentially  zero,  the 
rate  of  hydrogen  diffusion  through  the  oxide  layer  is  expected  to  be  (for  a  diffusion  process) 
proportional  to  -  Ami^tjij/Am, .  Substituting  the  right-hand  side  of  eq.  (80)  for  M^r  into 


this  expression  gives 


dW, 


sr  ■  -  -< 


AmH(t)  -Mho[1  -®‘U] 


(M> 


va 


r  -  G<Tw>&n  (82a) 


where  a  is  the  constant  of  proportionality,  or  what  might  be  termed  the  diffusion  constant. 
Here  ks  is  taken  from  eq.  (70),  using  the  values  of  Aj  and  Qs  given  in  eq.  (71).  To  convert 
eq.  (82a)  into  units  suitable  for  use  in  eq.  (12b),  it  must  be  multiplied  by  1 .06  x  10”3  to  con¬ 
vert  to  s!ugs/?ta-sec  and  by  the  area  of  the  fuel  element  to  obtain  G(TW)  in  slugs/sec,  i  .e. 


G(TW)  =  1-06  x  10-3  S^[AmH(t)] 


(82b) 


Equation  (82a)  must  be  integrated  to  obtain  Am^(t).  This  integration  gives  (see  Appendix  II) 

r  i/a 

Amj-|(t)  =  exp  (2a  J  t/k,' )  J  -i — Pi  -  exp  (2a  J t/k,f] 

^  1 0 


2a 


Xks 


i/8[ 1  +  exp  (-2a3  A  k$)] 


+  2{t^a  +  a/X  k/a  )[l  +  exp  (-2a8/X  ks)]  exp(-Xt  +  2a  ) 

”  « 

+  2/2tt'  exp  (-a8/X  k|)  ^erf(yTr  +  a//xkJP  )  -  erf  (a/  J  Xk$'  )]  ►  . 

. 


(83) 
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The  value  of  a  must  be  determined  from  experimental  data.  Using  sample  fuel  element  sped- 
mens  (provided  by  ARF  and  SRI)  that  had  been  exposed  to  oxidation  and  simulated  re-entry 
tests,  values  of  Am^t)  were  determined  by  chemical  analysis  (See  Appendix  I).  With  the  re** 
ported  values  of  M^0  and  these  measured  values,  eq.  (83)  was  used  to  determine  a.  Because 
t  was  very  large  for  the  oxidation  data  that  were  available,  Am|4  was  nearly  equal  to  Mjj0. 

Thus,  the  scatter  in  the  data  on  M^0  casts  considerable  doubt  on  the  values  of  a  so  determined. 
The  value  of  a,  based  upon  this  limited  data,  was  estimated  to  be 

a  =  0. 10  e”1  -84  x  1°VT  grams/cm3  -  min.  (84) 

Although  this  expression  may  not  be  very  accurate  in  general,  the  effects  involved  are  so  impor¬ 
tant  that  it  should  be  used  until  more  accurate  data  are  available.  It  would  not  be  difficult  to 
determine  an  accurate  value  of  a  by  resorting  to  experimental  methods,  e.g.,  weight  gain 
measurements  in  an  oxidizing  atmosphere. 

The  value  of  MHo 's  determined  by  dividing  the  fuel  element  area  by  the  mass  of  the 
total  hydrogen  contained  in  the  fuel  element. 

G.  DETERMINATION  OF  Hf,  Hv,  Twf,  Fc,  Fv,  HdaG,  AND  HC/  daG 

1 .  Heat  of  Fusion,  Hf,  and  Heat  of  Hydration,  HdaG 

These  quantities  can  be  taken  directly  from  the  report  by  Elliott.83  According  to  this 
report,  the  fuel  element  material  melts  at  3825°R,  and  the  heat  of  fusion  is 

Hf  »  90  Btu/lb  =  2.9  x  K^Btu/slug.  (85) 

Also,  the  report  gives  the  hydrogen  release  energy  as 

HdeG  =  -610Btu/,b  of  Zr'u’H  85  1-13  x  ICfBtu/slug.  (86) 
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2 .  Heat  of  Vaporization;  Hv  and  Twf 


Although  there  is  some  question  as  to  where  the  data  originated,  the  heat  of  vapori¬ 
zation,  according  to  a  NAAI  report,33  is 

Hv  =  180  Btu/lb  =  5.8  x  10*  Btu/slug.  (87) 

This  report  implies  that  the  vaporization  takes  place  at  a  temperature  of  4460°R,  hence,  the 
appropriate  value  of  Twf  would  be 

Twf  =  4460°F.  (88) 

A  calculation  would  then  have  to  be  made  to  see  if,  in  accordance  with  eq.  (11),  the  fuel 
element  temperature  would  remain  constant.  If  it  would,  then  the  fraction  of  the  material 
that  is  vaporized  is 

Fv  =  1.  <»> 

If,  under  these  conditions,  the  fuel  element  temperature  drops,  as  determined  by  eq.  (11), 
then  the  vapor  pressure  of  the  melted  material  will  have  to  be  known  unless  it  can  be  shown 
that  the  fuel  element  breaks  up  at  a  temperature  that  is  less  than  Twf.  All  presently  known 
data  indicate  that  the  fuel  element  material  will  remain  intact  until  the  metal  at  the  oxide 
layer-metal  interface  melts.  For  this  reason,  the  thermal  expansion  of  the  fuel  element  mate¬ 
rial  has  not  been  considered  in  the  foregoing  treatment. 

It  can  readily  be  shown  from  kinetic  theory  that  the  particle  flux,  leaving  a  surface 
whose  vapor  pressure  is  Pv  In  millimeters  of  Hg,  is  given  by  the  relation 

Flgx .  .II9Pv[£W*a]>iWft3^c  <«*> 

where  Mj«g  9ram  molecular  weight  of  the  evaporating  material,  in  the  case  of  zirconium. 


49 


91 .22.  Equation  (90)  can  be  multiplied  by  S  and  divided  by  g  and  then  substituted  Into 
the  left-hand  side  of  eq.  (11)  or  (13),  as  appropriate,  to  keep  track  of  the  transpiration 
cooling  effect  on  dTyy/dt  and  hence  on  Tw.  That  is. 


1  dW  0.119S 

~r~ 


,|U9,A*C  ■ 


(90b) 


Equation  (90)  can  also  be  employed  to  show  that  there  will  be  no  dearth  of  oxygen 
present  in  the  boundary  layer  for  oxidation  purposes  during  re-entry  simply  by  setting 
MdeG  =  32,  the  molecular  mass  of  oxygen,  and  using  the  local  boundary  layer  pressure  for 
Pv.  Such  a  calculation  shows  that  the  flux  available  (at  altitudes  below  300,000  ft)  is  cer¬ 
tainly  much  greater  than  can  be  used  in  the  oxidizing  process. 


3.  Heat  of  Combustion  or  Recombination,  Hc  jeQ ,  and  Fc 

There  is  a  possibility  that  the  evolved  hydrogen  atoms  will  recombine  and  then  bum 
with  the  oxygen  in  the  boundary  layer.  If  the  hydrogen  recombines,  the  energy  release  is 
approximately  3  x  10*  Btu/slug  (w  104  k  cal/ mole  of  Hg  ).  This  number  is  to  be  used  as 
Hc,  deG  for  the  hydrogen  recombination.  If  the  hydrogen,  after  recombination,  reacts  with 
the  oxygen  in  the  boundary  layer,  the  net  energy  release  will  be  approximately  3.5  x  1  (f 
Btu/slug  of  of  Hs  (120  k  cal/ mole  of  H* ) .  Since  the  activation  energy  for  the  formation  of 
Hs  from  H  atoms  is  zero,  it  follows  that  the  net  heat  energy  release  would  be  the  sum  of 
these  last  two  energies  (i.e.,  approx.  6.5  x  10# Btu/slug  of  hydrogen  used)  if  the  reaction 
went  straight  through  to  form  water  without  the  hydrogen  recombination  step. 

In  view  of  the  facts  that  the  burning  rate  of  hydrogen  in  oxygen  is  less  than  Mach  1 , 
the  boundary  layer  is  very  thin,  and  the  activation  energy  for  HsO  production  is  104  k  cal/ 
mole  of  HsO,  it  is  reasonable  to  assume  that  the  hydrogen  will  diffuse  well  into  the  boundary 
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layer  before  ary  appreciable  portion  of  it  goes  into  the  production  of  water.  Thu*,  the 
addition  of  the  combustion  energy  of  hydrogen  to  the  total  energy  content  of  the  boundary 
layer  is  expected  to  make  only  a  small  contribution  to  this  total  energy  content.  Conse¬ 
quently,  it  seems  reasonable  to  assume  that  Fc  can  be  considered  to  be  very  small  for  the 
combustion  of  Hs  and  that  the  Hacombustion  will  make  a  negligible  contribution  to  the  total 
heat  flux  into  the  fuel  rod  during  re-entry. 

The  case  for  neglecting  the  effects  due  to  the  heat  of  recombination  of  hydrogen  during 
the  re-entry  of  the  fuel  element  is  not  nearly  so  clear  cut  as  is  the  case  for  neglecting  the 
combustion  of  the  hydrogen.  None  of  the  data  available  to  the  authors  provides  information 
from  which  Fc  for  hydrogen  recombination  can  be  estimated.  In  view  of  the  fact  that  the  re¬ 
combination  energy  of  the  total  quantity  of  hydrogen  contained  in  the  fuel  element  is  about 
5y  500  Btu!s  (based  on  the  formula  ZrHl  8B-  U.0433)/  it  is  almost  as  large  as  the  total  aero¬ 
dynamic  heat  energy  available  to  the  fuel  rod  during  re-entry.  The  latter  quantity  is  only 
some  6,000  Btu's 

There  is  the  possibility  that  some  of  the  evolved  hydrogen  would  combine  at  grain 
boundaries  or  in  cracks  within  the  fuel  element,  and  some  of  it  might  combine  beneath  the 
oxide  layer.  Again,  some  of  it  might  be  evolved  before  the  stainless  steel  cladding  on  the 
fuel  element  has  burned  off  and,  thus,  recombine  between  the  fuel  material  and  the  cladding. 
However,  this  latter  possibility  seems  remote  because  the  stainless  steel  cladding  is  insulated 
from  the  fuel  material  by  a  ceramic  which  is  a  poor  heat  conductor. 

It  is  clear  from  the  foregoing  discussion  that  additional  research  is  required  before  an 
understanding  of  what  happens  to  the  hydrogen  in  the  fuel  element  can  be  realized. 
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IV.  SUMMARY  STATEMENT 


Although  the  effect  of  the  fuel  element  cladding  has  not  been  considered  here,  it  is 
obvious  that  the  same  treatment  is  applicable.  The  treatment  of  this  part  of  the  problem  that 
was  presented  by  NAAl®8  should  not  prove  grossly  in  error  when  used  with  the  foregoing  re-entry 
heating  equations. 

If  the  effects  due  to  the  recombination  and  combustion  of  the  hydrogen  content  of  the 
fuel  element  can  be  neglected  (it  seems  likely  that  the  effects  due  to  combustion, at  least,  can), 
then  bounds  of  the  upper  temperature  that  the  fuel  element  will  attain  during  re-entry  can  be 
calculated  w"th  a  high  degree  of  confidence  for  the  following  three  cases:  1)  no  transpiration 
cooling;  2}  all  hydrogen  transpires  at  the  phase  change  temperature  (approx.  2040^R);  3)  all 
hydrogen  transpires  during  the  peak  heating  portion  of  the  re-entry  trajectory. 

In  the  case  where  no  hydrogen  transpires,  a  hand  calculation,  using  the  NAAI^data, 
indicates  that  the  entire  fuel  element  will  be  melted  by  re-entry  heating.  If  the  hydrogen  is  all 
evolved  at  a  fuel  rod  temperature  of  2040°R  (phase  change  temperature),  the  estimated  peak 
tempercture  of  the  fuel  element  is,  based  on  hand  calculations,  approximately  4000°R.  If  it 
all  evolves  during  the  peak  heating  portion  of  the  re-entry  trajectory  (starting  at  about  280QPR), 
hand  calculations  indicate  that  the  peak  temperature  of  the  rod  will  be  approximately  3000°R. 
The  incident  heat  flux  used  in  these  calculations  was  taken  from  curves  presented  in  ref.  33  for 
re-entry  from  an  altitude  of  300,000  feet  at  24,200  ft/sec  and  an  initial  trajectory  angle  of 
-0.25  degrees.  The  continuous  transpiration  theory  derived  in  the  foregoing  analysis  does  not 
lend  itself  readily  to  a  hand  calculation  and  was  not  attempted  by  the  authors. 

Exclusive  of  the  behavior  of  evolved  hydrogen,  the  values  of  the  various  parameters  in 
the  re-entry  equations,  as  given  above,  are  quite  good.  Although  one  hesitates  to  place 
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accuracy  limits  on  the  entire  mathematical  model,  it  seems  unlikely  that  the  cumulative  error 
resulting  from  errors  in  the  values  of  the  parameters  given  above  would  exceed  10%  for  tem¬ 
peratures  up  to  the  melting  point  of  the  fuel  material . 

The  hydrogen  transpiration  data  presented  above  is  quite  susceptible  to  error  because 
of  the  poor  quality  of  the  initial  data  that  was  available  to  GTC .  This  is  not  meant  to  reflect 
poor  quality  work  on  the  part  of  other  researchers  who  produced  this  initial  data,  because  the 
experiments  and  tasks  they  performed  on  the  sample  materials  supplied  to  GTC  for  analysis  were 
not  designed  to  produce  the  type  of  data  that  GTC's  analysis  showed  to  be  required. 

At  the  present  time,  the;e  is  no  data  available,  to  the  authors'  knowledge,  from  which 
information  on  the  magnitudes  of  possible  effects  due  to  hydrogen  recombination  can  be  deduced. 
Also,  the  vapor  pressure  data  on  the  zirconium-uranium  alloy  could  not  be  uncovered.  Thus, 
the  value  of  the  parameters  Fc  and  Fv  are  undetermined  at  this  time. 

It  should  be  borne  in  miH  that  the  chemical  reactions  deduced  in  this  study  will  not 
necessarily  obey  the  same  laws  once  the  metal  becomes  molten  under  the  oxide  layer.  At  this 
point,  the  metal  cations  may  diffuse  to  the  surface  while  the  oxide  is  diffusing  to  the  molten 
material.  As  has  been  shown  by  Littman,  the  oxide  forms  a  solution  with  the  metal  when  the 
metal  is  in  a  molten  state.44  If  there  is  no  cation  diffusion  through  the  oxide,  then  the  con¬ 
tinued  oxidation  of  the  molten  metal  would  surely  disrupt  the  oxide  surface  layer  because  the 
density  ratio37 of  Zr  to  ZrOg  is  approximately  2.  If  the  oxide  layer  is  thus  broken  off,  the 
portion  of  the  rod  material  that  is  in  a  molten  state  will  blow  off  in  drop  or  droplet  form,  in 
which  case  the  theory  of  Opik4*  determines  the  maximum  size  of  the  particles.  The  portion  of 
the  rod  that  is  in  a  molten  state  at  any  time  can  be  calculated  from  the  system  of  re-entry 
equations  derived  above  and  eq .  (14) .  Whether  or  not  further  droplets  will  be  released  after 
this  initial  release,  should  such  an  initial  release  occur,  will  depend  on  the  vapor  pressure, 
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the  surface  tension,  the  viscosity,  etc .  of  the  molten  material .  It  seems  likely  that  the  zir¬ 
conium  does  diffuse  to  the  surface  to  relieve  the  internal  pressure  built  up  by  the  absorbed 
oxide,  because  Littman  found  that  when  cubes  of  the  fuel  element  material  have  an  oxide  layer 
on  their  surface,  they  retain  their  shape  when  inductively  heated  even  though  the  metal  inside 
is  in  a  molten  state.  If  this  is  the  case,  then  the  fuel  element  will  re-enter  intact  unless  break¬ 
away  oxidation  or  angular  acceleration  forces  due  to  tumbling  and  spinning  cause  the  oxide 
layer  to  break  off.  Although  this  report  shows  the  point  of  breakaway  oxidation  for  anion  diffu¬ 
sion,  this  value  is  not  expected  to  be  valid  when  both  anion  and  cation  diffusion  take  place 
simultaneously.  The  accelerating  forces  can  be  calculated  from  the  foregoing  equations  of 
angular  motion.  Should  it  turn  out  that  the  material  below  the  oxide  layer  becomes  molten, 
that  the  oxide  layer  is  thin,  and  that  the  angular  acceleration  forces  are  large,  the  oxide  layer 
could  be  broken  to  expose  the  molten  surface  to  viscous  drag  shearing  forces  which  will  produce 
molten  drops.  All  this  will  remain  speculative  until  the  actual  calculations  are  performed. 
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V.  REMAINING  RESEARCH  TO  BE  ACCOMPLISHED 


At  this  point  in  the  entire  SNAP  fuel  element  re-entry  program,  it  is  quite  clear 
what  remains  to  be  done  to  ascertain  whether  or  not  the  fuel  elements  will  burn  up  during 
re-entry.  This  could  be  accomplished  in  steps  as  follows: 

First,  the  research  necessary  to  determine  Fc  and  Fv,  together  with  a  better  value  of 
a  would  be  carried  out.  Then,  using  the  foregoing  analytical ,  mathematical  model,  a  machine 
calculation  would  be  performed  to  determine  if  any  substantial  portion  of  the  fuel  element  will 
reach  the  melting  point.  If  it  will  not,  then  the  fuel  element  will  not  burn  up  during  re-entry. 

It  is  apparent  that  the  upper  bound  heat  input  case  (transpiration  cooling  plus  hydrogen  recom¬ 
bination)  should  be  employed  in  the  first  calculation.  Thus,  if  this  calculation  were  to  show 
that  the  fuel  element  would  burn  up,  there  is  strong  evidence  that  it  could  be  made  to  burn  up 
d j'ing  /e-enny  . 

The  f  na!  stop  would  be  to  try  to  find  out  how  much  of  the  fuel  element  would  be  ablated 
awoy  Th'-s  would  require  a  detailed  knowledge  of  the  amount  of  material  that  blows  off  in 
droplet  form  and  from  which  part  of  the  fuel  rod  this  material  is  released  so  that  the  aerodynamic 
parameter  changes  could  be  accounted  for.  This  assumes  that  Fv  $  1 .  Should  research  show  that 
Fv  =  1 ,  then  the  foregoing  calculation  would  be  conclusive.  In  this  latter  case,  a  good  approxi¬ 
mation  to  the  degree  of  burn-up  could  be  obtained  by  assuming  that  the  material  evaporates  uni¬ 
formly  from  the  surface  of  the  rod  and  by  using  the  same  set  of  equations  given  above  with 
diminishing  diameter.  The  result  of  this  calculation  would  be  the  minimum  size  of  the  fuel  rod 
that  could  survive  the  re-entry.  If  this  is  zero,  then  there  is  little  to  worry  about  from  a  safety 
point  of  v:ew.  It  it  is  not  zero,  then  the  degree  of  hazard  associated  with  the  remaining  material 
could  be  specified  by  safety  personnel. 
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APPENDIX  I 


CHEMICAL  ANALYSIS 

This  appendix  contains  a  general  discussion  of  the  behavior  of  zirconium  and 
zirconium-uranium  alloys  in  various  gases.  The  chemical  analyses  and  results  obtained 
by  GTC  are  also  presented .  A  description  of  the  samples  of  fuel  element  materials  which 
GTC  received  during  the  contract  is  presented  in  tabular  form. 

A.  REVIEW  OF  PREVIOUS  STUDIES  ON  OXIDATION  OF  ZIRCONIUM  AND  THE  EFFECTS 
OF  O, ,  N„ ,  Ha  ON  ZIRCONIUM 

Since  preliminary  X-ray  data  indicated  that,  on  the  oxidation  specimens  provided  to 
GTC  by  the  Armour  Research  Foundation,  only  ZrOB  was  present  in  the  outer  layer  of  the 
specimen,  the  decision  was  made  to  survey  previous  investigations  on  the  kinetics  of  the  oxi¬ 
dation  ot  zirconium. 

The  reaction  of  zirconium  with  oxygen  at  various  temperatures  has  been  studied  by 
several  investigators.  A  comparison  of  the  results  shows  discrepancies  concerning  the  order 
of  the  rate  law  (cubic,  parabolic,  or  linear)  which  best  fits  the  oxidation  curves.  Hayes  and 
Roberson34  undertook  a  qualitative  study  of  this  subject  by  heating  panels  of  Zr  in  various 
gases.  They  found  that  above  800°C ,  oxygen  diffuses  into  the  grain  structure  of  the  metal 
when  heated  in  an  oxygen  atmosphere  at  1  atm  pressure.  It  appears  probable  that  the  oxygen 
penetration  into  the  grain  boundaries  takes  place  at  862°C,  when  the  close-packed  hexagonal 
alpha  phase  transforms  to  the  body  centered  cubic  beta  form.  The  two  percent  volume  change 
and  phase  transformation  certainly  provides  the  opportunity  for  intergranular  oxygen  diffusion. 
They  also  found  that  heating  to  900°C  in  moist  air  (percentage  of  water  vapor  not  given)  pro¬ 
duced  a  15  to  18  percent  volume  increase  of  the  specimen,  which  they  attributed  to  the  oxide 
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formation.  However,  heating  in  dry  air  or  oxygen  showed  no  perceptible  change  in  size  . 
Also,  heating  to  temperatures  above  1000°C  in  nitrogen  produced  a  nitride  film,  and  nitro¬ 
gen  diffused  into  the  metal,  but  the  nitride  film  formed  was  permeable  to  oxygen  at  1000°C . 

Hayes  and  Roberson  came  to  the  conclusion,  that  when  heated  in  air,  Zr  undergoes 
simultaneous  attack  by  both  oxygen  and  nitrogen.  The  nitride  was  found  to  be  unstable  in  the 
presence  of  oxygen,  and  decomposed  to  a  voluminous  oxide. 

Gulbransen  and  Andrew36  studied  the  reaction  on  zirconium  foil  specimens  between 
200°C  and  425°C .  They  reported  that  a  modified  parabolic  rate  law  fits  the  data.  The  tem¬ 
perature  dependence  of  the  oxidation  rate  as  a  function  of  film  thickness  showed  a  constant 
relationship.  This  fact  suggests  that  one  mechanism  is  controlling  the  rate  over  a  wide  thick¬ 
ness  range . 

In  another  study,  Gulbransen  and  Andrew36  found  that  between  400°  and  800°C,  thp 
method  of  surface  preparation  influenced  the  reaction  kinetics.  Mechanically  polished  foils 
obeyed  the  cubic  rate  law  and  reacted  faster  than  chemically  polished  foils,  which  obeyed 
the  parabolic  rate  law. 

A  later  study  by  Charles,  Barnartt,  and  Gulbransen,37 on  the  prolonged  oxidation  of 
zirconium  foils  up  to  times  of  500  hours  and  temperatures  of  350°  to  450°C,  confirmed  the 
cubic  rate  law  for  mechanically  polished  specimens.  However,  there  is  sufficient  reason  to 
suspect  that  the  size  and  shape  of  the  zirconium  samples  may  influence  the  reaction  kinetics. 
Belle  and  Mallett38  studied  the  oxidation  reaction  on  machined  cylindrical  specimens  4cm 
long  by  0.7cm,  and  5  cm  long  by  1  .4cm.  They  found  that  a  cubic  rote  law  fit  the  data 
taken  from  575°C  to  950°C  at  one  atmosphere  of  pure  oxygen,  while  on  foil  specimens  in 
approximately  the  same  temperature  range,  several  other  investigators  observed  a  parabolic 
rate  law  fit. 
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There  ts  also  evidence  that  the  particular  rate  law  which  the  reaction  follows  de¬ 
pends  on  the  temperature  range.  In  a  recent  study  by  Kofstad,39  zirconium  was  oxidized 
under  conditions  of  linearly  increasing  temperature.  It  was  shewn  that  between  6509  and 
950°C,  the  cubic  rote  law  was  obeyed,  and  between  950°  and  1100°C,  the  parabolic  rate 
law  fit  the  data.  This  is  in  close  agreement  with  the  analysis  of  the  ARF  data  presented  in 
the  text  of  this  report. 

Some  investigation  has  also  been  done  on  the  effect  of  pressure  on  the  oxidation  of 
zirconium.  Gulbransen  and  And  rev/,  among  others,  came  to  the  conclusion  that  pressure  has 
little  or  no  effect  on  the  reaction  of  zirconium  with  oxygen. 

The  reaction  of  zirconium  with  nitrogen  has  been  shown  to  proceed  at  a  much  slower 
rate  than  has  been  found  for  the  reaction  with  oxygen.  Several  studies  have  been  made  (see 
references  34  and  35),  and  most  investigators  agree  that  the  parabolic  rate  law  is  consistent 
with  the  data  at  the  higher  temperatures.  The  product  of  the  reaction  in  nitrogen  at  all  tem¬ 
peratures  has  been  identified  as  ZrN .  Gulbransen  and  Andrew34  noticed  that  small  amounts 
of  oxygen  in  the  nitrogen  noticeably  accelerate  the  rate  of  the  reaction.  Drounieks40  found 
that  at  temperatures  from  862°  to  1043°C,  the  nitrogen  diffuses  into  the  metal  as  well  as  pre¬ 
cipitates  a  nitride  film  on  the  surface. 

Some  work  has  also  been  done  on  the  reaction  of  zirconium  and  zirconium  alloys  with 
air.  In  general,  the  reaction  rate  is  more  rapid  in  air  than  in  either  oxygen  or  nitrogen  alone. 
In  this  case,  tho  solubilities  of  oxygen  and  nitrogen  in  the  metal  will  form  solid  solutions  of 
the  gases  in  the  metal  phase  as  well  as  phase  formation  of  ZrOfc  and  ZrN .  According  to 
Gulbransen,  for  the  zirconium  oxygen  reaction,  the  overall  rate  of  oxygen  absorption  still 
remains  parabolic,  if  the  rates  of  both  the  penetration  of  the  gas  Into  the  metal  and  the  growth 
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of  the  oxide  layer  are  diffusion  controlled,  and  the  Effusion  constants  remain  constant. 

The  same  type  of  reasoning  can  be  applied  to  the  nitrogen  reactions. 

Gulbransen  and  Andrew36  also  studied  the  rate  of  hydrogen  diffusion  into  zirconium 
at  300°C .  They  found  that  at  constant  te  nperature  the  reaction  rate  could  be  expressed  as 


dW 

3T  = 


where  p  is  the  pressure  of  the  hydrogen  atmosphere  and  A  is  a  constant.  According  to  theory, 
this  rate-pressure  relation  is  an  indication  of  activated  dirtusion  with  dissociation.  Therefore, 
the  correlation  shows  that  the  diffusion  of  hydrogen  into  zirconium  is  mainly  by  atoms  of  hydro¬ 
gen.  Additional  information  on  the  dissociation  rates  of  hydrogen  can  be  derived  from  the  data 
presented  by  Atkins..41  These  data  show  that  the  equilibrium  concentration  of  hydrogen  in 
Zr;H:U  is  very  nearly  independent  of  the  external  hydrogen  pressure,  varying  from  the  one- 
eighth  to  the  one -fourteenth  power  of  the  pressure. 

The  crystal  structure  of  ZrOg  has  been  studied  by  several  investigators.  Zirconium 
dioxide  has  been  shown  to  exist  in  three  crystalline  modifications  —  monoclinic,  cubic,  and 

tetragonal  .*  The  monoclinic  form  which  has  been  found  by  most  investigators  is  a  close 

o 

packed  arrangement  with  oxygen  spacings  of  2.6  to  2.7  A.  A  consideration  of  the  ionic  radii 
of  Zr+4 , 0.90  and  O'*,  1 .36 A,  shows  that  the  oxygen  ions  could  occupy  interstitial 
positions  in  the  ZrOg  lattice .  Thus,  the  oxidation  rate  is  expected  to  be  controlled  by  the 


rate  of  diffusion  through  the  oxide  layer. 


B.  RESULTS  OF  CHEMICAL  AND  PHYSICAL  ANALYSES 

The  problem  was  approached  by  using  various  techniques  of  chemical  and  physical 
analysis  to  obtain  data  from  certain  selected  ablation  and  oxidation  samples,  for  the  purpose 
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of  formalizing  and  substantiating  a  burn-up  model  of  the  zircontum-uronium-hydride  fuel 
element . 

The  samples  were  provided  to  GTC  by  Air  Force  contractors,  namely,  Vidya  Division 
of  the  Itek  Corporation,  Astropower,  Inc.,  Stanford  Research  Institute,  Battelle  Memorial 
Institute,  and  the  Armour  Research  Foundation.  As  the  samples  were  received  by  GTC,  they 
were  tabulated  along  with  the  simulated  conditions  under  which  they  were  tested  (see  Table  5). 
All  samples  which  had  been  prepared  for  metallographic  study  were  examined  under  a  microscope 
for  the  presence  of  cracks,  voids,  oxidation  layers,  disrupted  grain  boundaries  and  grain  struc¬ 
ture. 

The  samples  to  be  subfected  to  X-ray  and  chemical  analysis  were  chosen  on  the  basis 
of  the  test  conditions  to  which  they  were  exposed.  A  selection  was  made  of  twenty  samples 
to  be  chemically  analyzed,  and  five  samples  were  used  for  X-ray  analysis. 

For  the  determination  of  compounds  produced  by  the  effects  of  ablation  and  oxidizing 
conditions,,  a  determination  of  crystal  structures  was  made  by  X-ray  diffraction .  Although 
this  method  is  limited  to  qualitative  results,  the  identity  of  the  ma{or  compounds  which  are  in¬ 
volved  in  the  process  was  established. 

The  X-ray  data  were  taken  on  a  Siemens  X-ray  diffraction  unit,  using  a  30  kilovolt 
source  with  20  ma  K-a  Cu  radiation .  The  patterns  were  scanned  at  0.5°per  minute  from  15° 
to  110°  with  a  2%  statistical  error.  The  peaks  were  checked  against  the  1962  NBS  card  data 
f'le  for  the  following  compounds: 

1 .  Both  aand  S  phase  zirconium  dioxide  (ZrOg) 

2  Uranium  metal  (U) 

3  .  Uranium  nitride  (UN) 

4  fi  Uranium  metal  (8U) 

5 .  Di uranium  pentaoxide  (Us Ob) 

6  .  Three  phases  of  triuranium  octooxide  (UjC^) 
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7.  a  and  0  uranium  dioxide  (UOg) 

8 .  Uranium  trioxide  ( U03  ) 

9.  Uranium-zirconium  (U-Zn) 

10.  Zirconium  nitride  (ZrN) 

The  results  of  this  analysis  are  given  in  Table  3. 


TABLE  3 

X-RAY  DIFFRACTION  RESULTS 


GTC  Nr. 

Material 

Test  Temp. 
(°C) 

Air-Argon 
Atmosphere 
(%  Air) 

Sample  Condition 

Compounds  Identified 

26 

Zr-U-H 

1100 

50 

Powdered  outer  layer 

ZrCfe,  possibly  uranium 
metal 

20 

Zr-U-H 

800 

50 

Surface  film 

Results  cannot  be  identi¬ 
fied,  possibly  because 
of  interference  due  to 
substrate  layer 

25 

Zr-U 

1100 

50 

Surface  film 

Results  cannot  be  positively 
identified,  but  ZrOg 
indicated 

28 

Zr-U-H 

1100 

100 

Surface  film 

ZrOB  only 

32 

Zr-U-H 

1400 

50 

Powdered  outer  layer 

Very  strongly  ZrQj 
(24  peaks  matched) 

The  chemical  analysis  of  the  samples  was  made  by  a  vacuum  fusion  technique.  This  type 
of  analysis  gives  the  weight  percentages  of  oxygen,  nitrogen,  and  hydrogen  remaining  in  the 
samples  after  testing  under  simulated  ablation  conditions.  Because  the  samples  are  all  inorganic 
substances  with  high  melting  points  and  a  quantitative  measure  of  the  combined  and  absorbed 
gases  is  necessary  to  supply  data  for  the  reaction  constants  and  hydrogen  evolution  rate,  the  high 
temperature  vacuum  technique  is  about  the  only  method  by  which  the  small  amounts  of  samples 
can  be  analyzed . 
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Where  possible,  certain  samples  were  analyzed  In  two  parts.  One  part  consisted  of 
the  surface  layer,  composed  mainly  of  ZrOa  and  the  other  part  was  the  center  core  of  the  sam¬ 
ple  minus  the  surface  layer.  This  technique  provided  information  on  the  amount  of  oxygen  and 
nitrogen  which  had  diffused  into  the  center  section  of  the  fuel  element. 

The  results  of  the  vacuum  fusion  analysis  for  the  twenty  samples  appear  in  Table  4. 

Since  S.R.I .  ran  samples  under  constantly  changing  conditions  comparable  to  re-entry 
trajectories,  only  generalizations  can  be  made  about  the  results.  It  appears  that  up  to  3000^R,a 
significant  amount  of  hydrogen  is  retained  by  the  sample.  It  also  is  apparent  that  the  build-up 
of  an  oxidation  layer  on  the  surface  presents  a  definite  barrier  to  the  evolution  of  hydrogen, 
even  up  to  3500°R.  Furthermore,  it  is  apparent  from  l.l.T.R.I.  samples  that  the  hydrogen  is 
retained  in  the  samples  up  to  1 100°C  regardless  of  the  amount  of  oxide  formation  on  the  surface 
of  the  sample.  In  most  cases,,  at  temperatures  of  1100°C  and  upward,  a  large  amount  of  oxygen 
is  present  in  the  surface  layer,  and  some  of  the  oxygen  has  managed  to  diffuse  into  the  center 
core  of  the  sample.  Judging  from  the  increases  in  the  nitrogen  content  of  some  of  the  samples 
nitrogen  combination  or  absorption  plays  some  role  in  the  chemical  kinetics  of  the  ablation  pro¬ 
cess,  although  its  role  is  not  nearly  as  significant  as  that  of  oxidation  reactions.  Unfortunately, 
no  conclusions  can  be  derived  concerning  the  effects  or  probability  of  a  hydrogen -oxygen  re¬ 
action  from  the  data  presented .  In  certain  cases,  the  loss  of  hydrogen  may  be  responsible  for  a 
decrease  in  oxidation,  while  in  other  cases  it  may  promote  oxidation. 

An  interpretation  of  the  analytical  results  in  terms  of  specific  reactions  is  difficult  be¬ 
cause  the  effects  are  of  an  integrated  nature  and  therefore  difficult  to  separate  experimentally. 
However,  the  data  were  most  helpful  in  determining  reaction  rate  constants,  in  establishing  the 
order  of  the  oxidation  rate  law,  and  in  determining  the  diffusion  coefficient  for  the  hydrogen  loss. 
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TABLE  4 


VACUUM  FUSION  ANALYSIS  OF  SAMPLES 


Gas  Analysis  (Wt.%) 

GTC  No.  Material  Test  Temp ■  Atmosphere  Sample  Condition  Oxygen  Nitrogen  Hydrogen 

Blank  (A.I.  Material) 


0  (blank) 

Zr-U-H 

None 

None 

Chunk 

.064 

.019 

1 .64 

Stanford  Research  Samples 

51 -H 

Zr-U-H 

3500°R 

Air 

Pyramid 

.139 

.028 

.073 

53-H 

Zr-U-H 

3840°R 

Air 

Pyramid 

.420 

.169 

.020 

56-H 

Zr-U 

3590°R 

Air 

Half-Pyramid 

2.800 

.800 

.003 

56-H 

Zr-U-H 

3550°R 

Air 

Half-Pyramid 

3.980 

070 

.148 

40-H 

Z--U-H 

2940°R 

Air 

Pyramid 

.065 

.012 

476 

l.l.T.B 

' .  1 .  (Formerly  Armour  Research)  Samples 

12 

Zr-U-H 

500°C 

Ar-10%  Air 

1/2  Disk 

.289 

.017 

1.078 

16 

Zr-U-H 

500°C 

100%  Air 

1/2  Disk 

.386 

.109 

1.000 

18 

Zr-U-H 

800°C 

Ar-10%  Air 

1/2  Disk 

1.022 

.067 

.930 

21 

Zr-U 

800°C 

100%  Air 

Fragments 

.146 

.002 

.005 

22 

Zr-U-H 

800°C 

100%  Air 

Fragments 

.231 

.002 

1.130 

24(a) 

Zr-U  SI 

noo°c 

Ar-10%  Air 

1/2  Center  core 

1.050 

.062 

.008 

24(b) 

Zr-U-H 

1100°C 

Ar-10%  Air 

Surface  layer 

4.028 

.506 

.014 

27(a) 

Zr-U 

noo°c 

100%  Air 

Fragments  of 

center  core 

.048 

.009 

.002 

27(b) 

Zr-U 

1100°C 

100%  Air 

Surface  layer 

.638 

.613 

.014 

28 

Zr-U-H 

11C0°C 

100%  Air 

Fragments 

.613 

.025 

.995 

30(a) 

Zr-U-H 

1400°C 

Ar-10%  Air 

Core 

3.160 

.127 

.027 

30(b) 

Zr-U-H 

1400°C 

Ar-10%  Air 

Surface  layer 

11.992 

2.202 

.015 

33(a) 

Zr-U 

I400°C 

100%  Air 

Core  fragments 

1.350 

.025 

.008 

33(b) 

Zr-U 

1400°C 

100%  Air 

Surface  layer 

8.930 

.006 

.010 

Vidya  Sample 

5 

Zr-U-H 

Unknown 

Air 

Surface  layer  from 

3.695 

.354 

.035 

exposed  surface 
of  cylinder 
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TABLE  5 

TABULATION  OF  SAMPLES 
A.  Vidyo  Samples 


All  models  were  one  inch  diameter  cylinders  of  hydrided  fuel  rod  material. 
Atmosphere;  dry  air. 


Sample 

G.T.C. 

i  No. 

Vidya 

Avg.  cold  Stagnation 
wall  heat  enthalpy 
rate® .  Hst 

(BTU/fr-sec)  (BTU/lb) 

Stagnation 
pressure  at 
model  surface 
Pst  (psia) 

Model 

weight 

change 

(grams) 

Run 

duration 

(seconds) 

1 

538 

104 

6380 

14.7 

0.000 

60 

la 

544b 

845 

3770 

16.2 

-64.49 

59.5 

2 

539 

605 

2280 

19.0 

-.350 

60 

3 

550 

868 

3180 

19.9 

-59.40 

60 

4 

55! 

163 

4930 

14.7 

+.016 

121 

5 

552 

582 

3980 

16.0 

+.262 

62.3 

6 

553 

93 

7080 

14.7 

+.020 

121.3 

7 

554 

574 

3810 

16.0 

-49. 187 

92 

8 

555 

453 

4680 

15.5 

-.062 

61 

9 

556 

876 

3230 

19.0 

+.341 

70.4 

10 

557 

435 

4630 

15.5 

+.228 

120.2 

Notes:  a.  Heat  rate  averaged  over  the  surface  of  a  flat  face  or  hemispherical  surface 
of  a  cooled  copper  calorimeter  of  the  same  size  and  geometry  as  the  test 
model. 

b.  Rerun  of  run  538  with  nitrogen  as  test  fluid. 
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TABLE  5  (Continued) 

B.  Illinois  Institute  of  Technology  Research  Institute  Samples 


In  each  case,  the  time  of  run  was  3  hours,  at  1  atmosphere  pressure. 


G.T.C.  Sample 
No. 

Composition 

Air-Argon 
Atmosphere  (%  Air) 

Heat  flux 
or  temp.(°C) 

Comments 

11 

Zr-U 

10 

500 

Grey  outside,  thin  layer 

12 

ZrH-U 

10 

500 

1  mil  thick  outside  layer 

i3 

Zr-U 

50 

500 

Grey  outside  layer 

14 

ZrH-U 

50 

500 

Small  cracks  around  edges 

15 

Zr-U 

100 

500 

Needle-like  grain  structure 
but  no  oxide  formation 

16 

ZrH-U 

100 

500 

Shows  distorted  grain  boundaries 

17 

Zr-U 

10 

800 

Grey-white  outside  layer 

18 

ZrH-U 

10 

800 

Definite  oxide  layer,  distorted 
grain  boundary  near  edges 

19 

Zr-U 

50 

800 

Grey-white  outside  layer 

20 

ZrH-U 

50 

800 

1-1/2  mil  thick  outside  Icyer, 
disruption  of  grain  boundaries 

21 

Zr-U 

100 

800 

Grey-white  outside  layer 

22 

ZrH-U 

100 

800 

Random  grain  structure  near 
edges  showing  surface  reaction 

23 

Zr-U 

10 

1100 

Oxide  layer  1/32  inch  deep 

24 

ZrH-U 

10 

1100 

Material  flakes  off  surface 

25 

Zr-U 

50 

1100 

Thin  white  oxide  coat,  grey 
beneath 

26 

ZrH-U 

50 

1100 

Silver-grey  oxide  layer,  gold 
flecked  beneath 
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TAP!  E  5  (Continued) 

B.  illinoii  Institute  of  Technology  Research  Institute  Samples  (Continued) 


G.T.C.  Sample 
No. 

Composition 

Air-Argon 
Atmosphere  (%  Air) 

Heat  flux 
or  temp.  (°C) 

Comments 

27 

Zr-U 

100 

1100 

White  oxide  layer  on  surface 

2d 

ZrH-U 

100 

1100 

Black  pitted  surface,  cut  and 
pol  ished  on  one  side 

29 

Zr-U 

10 

1400 

Heavy  oxide  coating 

30 

ZrH-U 

10 

1400 

Heavy  oxide  coating 

31 

Zr-U 

50 

1400 

Heavy  oxide  coating 
(powdery) 

32 

ZrH-U 

50 

1400 

Oxidized  all  the  way  through 

33 

Zr-U 

100 

1400 

Oxidized  all  the  way  through 

34 

ZrH-U 

100 

1400 

Cracked  and  spalling 
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TABLE  5  (Continued) 

C.  Stanford  Reseorch  Institute  Somples 
"""“"Air  Atmosphere 


. Meat  ^jux 


Sample  No. 

G. T. C.  S.  R.  1 . 

Trajectory 

Comp. 

Av0-  o 

(BTU/fr 

Actual 

-sec) 

Temp. 

(°R) 

System 
Pressure 
(in.  of  Hg.) 

Time 

(sec.) 

Comments 

37  a-1 

alloy 

— 

— 

— 

— 

— 

Not  tested 

38  la-2 

alloy 

110 

165 

2350 

0.1 

100 

Grey-white  layer. 

150 

225 

2560 

0.2 

100 

cracked  on  one  c>dge . 

265 

395 

2960 

0.6 

100 

365 

545 

3140 

2.0 

50 

215 

120 

2780 

4.0 

50 

45 

25 

2240 

2.6 

50 

39  la-4 

hydride 

110 

165 

2290 

0.1 

100 

Grey  coating,  eroded 

150 

225 

2380 

0.2 

100 

on  edges. 

265 

395 

2830 

0.6 

100 

365 

545 

3120 

2.0 

50 

215 

120 

2890 

4.0 

50 

45 

25 

2260 

2.6 

50 

40  la-9 

hydride 

110 

165 

2260 

0.1 

100 

Grey  coating,  eroded 

150 

225 

2350 

0.2 

100 

on  edges. 

265 

395 

2740 

0.6 

100 

365 

545 

2940 

2.0 

50 

215 

120 

2600 

4.0 

50 

45 

25 

2220 

2.6 

50 

41  2a~6 

alloy 

365 

545 

3100 

0.8 

50 

Grey-white  coating. 

360 

540 

2995 

1.0 

50 

no  erosion. 

265 

395 

2650 

1.9 

50 

100 

150 

2260 

3.7 

50 

10 

15 

2095 

1.6 

50 

42  2a-9 

hydride 

365 

545 

2840 

0.8 

50 

Grey  coating,  edges 

360 

540 

2350 

1.0 

50 

eroded. 

265 

395 

2610 

1.9 

50 

100 

150 

2260 

3.7 

50 

10 

15 

2130 

1.6 

50 

43  3a- 1 

alloy 

700 

900 

3480 

4.5 

20 

Grey-whi  te  surface 

340 

510 

3360 

4.0 

20 

layer. 

160 

240 

2470 

4.4 

20 

40 

60 

2400 

3.9 

20 
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TABLE  5  (Continued) 

C.  Stanford  Reseorch  Institute  Samples 
"*  Air  Atmosphere  ((Hontinuea) 


Sample 

G.T.C. 

No. 

S.  R.  1. 
Trajectory 

Comp. 

Heat  Flux 
TTt 

Avg.  Actual 
(BTU/ft2 -sec) 

Temp. 

W 

System  Time 

Pressure  (sec.) 

( in.  of  Hg.) 

Comments 

44 

2b- 1 

alloy 

245 

360 

2565 

1.2 

50 

All  white  oxide 

450 

675 

3280 

8.8 

50 

covered,  ends  and 

720 

900 

3515 

36.8 

50 

edges  split. 

400 

600 

3190 

62.0 

50 

45 

3b- 1 

alloy 

250 

375 

2470 

10.4 

35 

Greyish  -wh  i  te ,  ends 

400 

600 

3300 

44.0 

35 

and  edges  split  open. 

240 

360 

2470 

50.4 

35 

46a 

a -4 

alloy 

45 

68 

2200 

* 

300 

Sectioned  and  potted. 

90 

145 

2350 

.04 

100 

shows  sma  1 1  amount  of 

200 

300 

2710 

.24 

100 

surface  reaction. 

340 

510 

3090 

1.40 

100 

140 

210 

2560 

2.40 

100 

46h 

a-o 

hydride 

45 

68 

2190 

* 

300 

Sectioned  and  potted. 

90 

145 

2280 

.04 

100 

shows  smal  1  amount  of 

200 

300 

2440 

.24 

100 

surface  reaction. 

340 

510 

2930 

1.40 

100 

140 

210 

2490 

2.40 

100 

47a 

a -7 

alloy 

45 

68 

2330 

* 

300 

Sectioned  and  potted. 

90 

145 

2090 

.04 

100 

shows  layer  of  surface 

200 

300 

3190 

.24 

100 

reaction. 

340 

510 

3500 

1.40 

100 

140 

210 

2800 

2.40 

too 

47h 

a -9 

hydride 

45 

68 

2200 

* 

300 

Sectioned  and  potted. 

90 

145 

2240 

.04 

100 

shows  layer  of  surface 

200 

300 

2510 

.24 

100 

reaction. 

340 

510 

3010 

1.40 

100 

140 

210 

2320 

2.40 

100 

48a 

la-3 

alloy 

110 

165 

2550 

0.1 

100 

Sectioned  and  potted. 

150 

225 

2670 

0.2 

100 

grain  boundaries dis- 

205 

395 

2980 

0.6 

100 

rupted  near  outer 

305 

545 

3120 

2.0 

50 

surface. 

215 

120 

2750 

4.0 

50 

45 

25 

2210 

2.6 

50 

*no  data 
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TABLE  5  (Continued) 

C.  Stanford  Research  Institute  Samples 
Air  Atmosphere  (Continued) 


Heat  Flux 

#  e 


Sample  No. 
G.T.C.  S.  R.l. 

Trajectory 

Comp. 

Avg.  Actual 
(BTU/fr  -sec) 

Temp. 

(°R) 

System  Time 

Pressure  (sec.) 
(in.  ofHg.) 

Comments 

48h 

la-l 

hydride 

110 

165 

2330 

0.1 

100 

Sectioned  and  potted , 

150 

225 

2350 

0.2 

100 

grain  boundaries  dis¬ 

205 

395 

2820 

0.6 

100 

rupted  near  outer 

305 

545 

3120 

2.0 

50 

surface. 

215 

120 

2890 

4.0 

50 

45 

25 

2260 

2.6 

50 

49a 

2a -7 

alloy 

365 

545 

2970 

0.8 

50 

Sectioned  and  petted. 

360 

540 

3010 

1.0 

50 

265 

395 

2670 

1.9 

50 

100 

150 

2240 

3.7 

50 

10 

15 

2095 

1.6 

50 

49h 

2a -8 

hydride 

365 

545 

2830 

0.8 

50 

Hydride  is  cracked 

360 

540 

2850 

1.0 

50 

along  back  edge. 

265 

395 

2615 

1.9 

50 

Sectioned  and  potted. 

100 

150 

2260 

3.7 

50 

10 

15 

2132 

1.6 

50 

50a 

3a -2 

alloy 

700 

900 

3500 

4.5 

20 

Sectioned  and  potted . 

340 

510 

3320 

4.0 

20 

160 

240 

2545 

4.4 

20 

40 

60 

2420 

3.9 

20 

50h 

3a -3 

hydride 

700 

900 

3140 

4.5 

20 

Hydride  cracked 

340 

510 

3010 

4.0 

20 

along  back  edge. 

160 

240 

2400 

4.4 

20 

40 

60 

2310 

3.9 

20 

51a 

4a -5 

alloy 

480 

720 

3715 

20.8 

20 

Sectioned  and  potted , 

140 

210 

2330 

4.4 

20 

alloy  has  surface  re¬ 

48 

72 

2110 

2.4 

20 

action. 

51h 

4a  -8 

hydride 

480 

720 

3500 

20.8 

20 

Sectioned  and  potted. 

140 

210 

2185 

4.4 

20 

hydride  is  cracked. 

48 

72 

2110 

2.4 

20 
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TABLE  5  (Continued) 

C.  Stanford  Rese arch  Institute  Samp  lei 
Air  Atmosphere  ((Continued) 


Heat  Flux 
•  •  _ 


Sample  No. 
G.T.C.  S.R.I. 

Trajectory 

Comp. 

Avg.  .Actual 
(BTU/ft2  -sec) 

Temp. 

(°R) 

System  Time 

Pressure  (sec.) 
(in. of  Hg) 

Comments 

52a 

b-3 

alloy 

60 

90 

2470 

0.2 

400 

Sectioned  and  potted, 

130 

195 

2600 

0.4 

200 

cracked. 

320 

480 

3000 

2.4 

200 

640 

900 

2650 

18.4 

100 

320 

480 

3020 

66.0 

100 

52h 

b-4 

hydride 

60 

90 

2470 

0.2 

400 

Sectioned  and  potted, 

130 

195 

2550 

0.4 

200 

cracked. 

320 

480 

2980 

2.4 

200 

640 

900 

2650 

18.4 

100 

320 

480 

2990 

66.0 

100 

53a 

b-6 

alloy 

60 

90 

2510 

0.2 

400 

Sectioned  and  potted, 

130 

195 

2650 

0.4 

200 

reaction  layer  around 

320 

480 

3155 

2.4 

200 

edges  of  specimen. 

640 

900 

3695 

18.4 

100 

320 

480 

3155 

66.0 

100 

53h 

b-8 

hydride 

60 

90 

2470 

0.2 

400 

Sectioned  and  potted, 

130 

195 

2615 

0.4 

200 

reaction  layer  around 

320 

480 

3245 

2.4 

200 

edges  of  specimen. 

640 

900 

3840 

18.4 

100 

320 

480 

3140 

66.0 

100 

54a 

lb-3 

alloy 

65 

100 

2310 

0.2 

300 

Sectioned  and  potted, 

215 

325 

2490 

1.2 

100 

definite  layer  around 

415 

625 

3120 

4.8 

100 

edges  of  sample. 

670 

900 

3640 

30.0 

100 

240 

360 

2565 

30.0 

100 

54h 

lb-4 

hydride 

65 

100 

2350 

0.2 

300 

Sectioned  and  potted, 

215 

325 

2435 

1.2 

100 

definite  layer  around 

415 

625 

3210 

4.8 

100 

edges  of  sample. 

670 

900 

3720 

30.0 

100 

240 

360 

2760 

30.0 

100 
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TABLE  5  (Continued) 

C.  Stanford  Research  Institute  Samples 
Air  Atmosphere  (Continued) 


#Heat  F|ux 


Sample  No. 

G.  T.  C.  S.  R.  1 . 

Trajectory 

Comp. 

*7  ^ 

Avg.  Actual 
(BTU/ft2  -»ec) 

Temp. 

(°R) 

System  Time 

Pressure  (sec.) 
(in. of  Hg.) 

Comments 

55a 

1b-5 

alloy 

65 

100 

2110 

0.2 

300 

Sectioned  and  potted, 

215 

325 

2760 

1.2 

100 

[ 

415 

625 

3660 

4.8 

100 

1 

670 

900 

3800 

30.0 

100 

240 

360 

2740 

30.0 

100 

55h 

lb-4 

hydride 

65 

100 

2145 

0.2 

300 

Sectioned  and  pitted 

i 

215 

325 

2185 

1.2 

100 

t 

415 

625 

3280 

4.8 

100 

670 

900 

3680 

30.0 

100 

l 

240 

360 

2650 

30.0 

100 

56a 

2b -2 

alloy 

245 

360 

2510 

1.2 

50 

Sectioned  and  potted 

450 

675 

3320 

8.8 

40 

720 

900 

3590 

36.8 

50 

1 

400 

600 

3210 

62.0 

50 

56h 

2b -4 

hydride 

245 

360 

2290 

1.2 

50 

Sectioned  and  potted 

450 

675 

2940 

8.8 

40 

hydride  cracked,  re¬ 

| 

720 

900 

3550 

36.8 

50 

action  layer  around 

f 

400 

600 

3110 

62.0 

50 

edges. 

1 

!  57a 

3b -2 

alloy 

250 

375 

2580 

10.4 

35 

Sectioned  and  potted 

400 

600 

3190 

44.0 

35 

slight  recovery  zone 

240 

360 

2435 

50.4 

35 

on  surface. 

57h 

3b -4 

hydride 

250 

375 

2330 

10.4 

35 

Sectioned  and  potted 

400 

600 

3010 

44.0 

35 

slight  recovery  zone 

240 

360 

2330 

50.4 

35 

on  surface. 

\ 
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TABLE  5  (Continued) 

D.  Armour  Reseorch  Foundotlon 
Tested  at  I  Atm.  Pressure 


Sample  No. 
G.T.C. 

Composition 

Type  of  test 
performed 

Atmos¬ 

phere 

Temp. 

<°R> 

Comments 

58 

ZrH-U 

Loss  of  fuel 
rod  integrity 

Argon 

3350 

5  pieces  showing  zones  of 
reaction. 

59 

ZrH-U 

Linear  thermal  Hydrogen 
expansion 

2790 

Long,  rod-shaped  piece. 

60 

Zr-U 

Linear  thermal  Argon 
expansion 

2790 

Long,  rod-shaped  piece. 

61 

ZrH-U 

Thermo  1 
conductivity 

Hydrogen  2710 

1"  rods,  material  flakes 
off. 

62 

Zr-U 

Thermal 

conductivity 

Argon 

2750 

1"  rods,  material  flakes 
off. 

63 

ZrH-U 

Emissivity 

Air 

2500 

Round  disks  about  1.343" 

64 

ZrH-U 

Emissivity 

Hydrogen 

2800 

diameter,  material  flakes 

65 

ZrH-U 

Emissivity 

Argon 

2800 

off. 

66 

Zr-U 

Specific  heat 

Argon 

2770 

Rod-shaped  pieces,  sur- 

67 

2  samples  ZrH-U 

Specific  heat 

Hydrogen 

2700 

face  reaction. 

68 

Zr-U 

Heat  of  fusion 

Argon 

3430 

Encapsulated  in  graphite. 

69 

ZrH-U 

Tensile  test 

Air 

Room 

temp. 

Broken  at  gauge  mark. 

70 

ZrH-U 

Tensile  test 

Air 

Room 

temp. 

Broken  in  two  places. 

71 

Zr-U 

Tensile  test 

Air 

Room 

temp. 

Broken  inside  gauge  mark. 

72 

Zr-U 

Tensile  test 

Argon 

2600 

Not  broken  but  partially 
melted  at  one  spot. 
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TABLE  5  (Continued) 

E.  Battalia  Instituta  Samples 


samples  subjected  to  rocket  exhaust 


Sample  No.  Comj^ssitlon  test  In  O2/  H2,  and  H2O  atmosphere  at  Comments 
G.T.  C.  3500°R.and  0.25  atm. pressure  unless 

otherwise  noted 


73 

ZrH-lOU  (A.  1.,  material) 

Melted  on  face  of  cylinder. 

74 

1 -ZrH-l.OU  -2Ba 

Sectioned  and  potted. 

75 

2-ZrH-10U-2Nb 

Sectioned  and  potted. 

76 

3-ZrH-10U-5Sn 

Sectioned  and  potted. 

77 

4-ZrH-10U-2Cu 

Sectioned  and  potted. 

78 

5-ZrH-10U-2V 

Sectioned  and  potted. 

79 

6-ZrH-10U-2Ca 

Sectioned  and  potted. 

80 

7-ZrH-10U-2Nb-2Ba 

Sectioned  and  potted. 

81 

8-ZrH-10U-2Sn-2Nb 

Sectioned  and  potted. 

82 

9-ZrH-10U-5Cr 

Sectioned  and  potted. 

83 

10-ZrH-10U 

Sectioned  and  potted. 

84 

1 1-ZrH-10U-Ba 

Sectioned  and  potted. 

85 

l2-ZrH-10U-2W 

Sectioned  and  potted. 

86 

13-ZrH-10U-2Ta 

Sectioned  and  potted. 

87 

14-ZrH-10U-2AI 

Sectioned  and  potted. 

88 

15-ZrH-10U-2Bi 

Sectioned  and  potted. 

89 

16-ZrH-10U-5Sb 

Sectioned  and  potted. 

90 

17-ZrH-10U  (A.I  .  material) 

Exposed  to  rocket  5  seconds, 
sectioned  and  potted. 

91 

18-ZrH2-U02 

Exposed  to  rocket  5  seconds, 
sectioned  and  potted. 

92 

19-ZrH2-UC 

Exposed  to  rocket  5  seconds, 
sectioned  and  potted. 

93 

20-ZrH2-UN 

Exposed  to  rocket  5  seconds, 
sectioned  and  potted. 

94 

2l-ZrH-10U-2Sr 

Exposed  to  rocket  5  seconds, 
sectioned  and  potted. 

95 

22 -ZrH-lOU -5Nb 

Exposed  to  rocket  5  seconds, 
sectioned  and  potted. 

96 

23-ZrH-10U-5Th 

Exposed  to  rocket  5  seconds, 
sectioned  and  potted. 

97 

24-ZrH-10U-2Ni 

Exposed  to  rocket  5  seconds, 
sectioned  and  potted. 

98 

25- ZrH-lOU -2Si 

Exposed  to  rocket  5  seconds, 
sectioned  and  potted. 

99 

26-ZrH-10U-2Cr 

Exposed  to  rocket  5  seconds, 
sectioned  and  potted. 

100 

27 -ZrH-lOU  (A.  1.  material) 

Exposed  to  rocket  2.5  seconds, 
sectioned  and  potted. 

101 

ZrH-10U-5Sn 

Not  tested. 
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TABLE  5  (Continued) 

E.  Bottelle  Institute  Samples  (Continued) 

~  [  All  samples  subjected  to  rocket  exhaust 

Sample  No.  tost  in  H,,  and  H2O  atmosphere  at 

G.T.C.  Composition  3500°R.  and  0.25  atm.  pressure  unless  Comments 

otherwise  noted 


102 

ZrH-10U-2AI 

Exposed  to  rocket  for  8  seconds, 
cylinder  ablated  on  one  end. 

103 

ZrH-10U-2W 

Exposed  to  rocket  for  8  seconds, 
cylinder  ablated  on  one  end  half¬ 
way  back. 

104 

ZrH-10U-2V 

Exposed  to  rocket  for  8  seconds, 
cylinder  ablated  on  one  end. 

105 

ZrH-10U-2Si 

Not  tested. 

106 

Unidentified  alloy  section 

Specimen  found  down  stream  after 
test. 

107 

ZrH-UN 

Section  of  isostatically  compacted 
specimens  found  down  stream. 

108 

ZrH-UC 

Section  of  isostatically  compacted 
specimens  found  down  stream. 

109 

Probably  ZrH-U02 

Section  of  isostatically  compacted 
specimens  found  down  stream. 

110 

Probably  ZrH-U02 

Section  of  isostatically  compacted 
specimens  found  down  stream. 

111 

Suspected  to  be  ZrH-UN 

Section  of  isostatically  compacted 
specimens  found  down  stream. 

112 

Suspected  to  be  ZrH-UN 

Section  of  isostatically  compacted 
specimens  found  down  stream. 

TABLE  5  (Continued) 

F.  Astropower,  Inc.  Samples 
Ail  Tatting  Performed  In  Arc  Plosmo  Facility 


Sanipl«  NoT  Partial  Pressures  of  Plasma  Gases 

G.T.C.  Astropower  Composition  Argon(qtfn.)  Oxypenfatm.)  Nitrogenfetm.)  Comments 


113 

64 

Zr-U 

1.000 

0 

0 

Mounted  particles. 

1 14 

65 

Zr-U 

.667 

0 

.333 

Mounted  particles. 

1 15 

66 

Zr-U 

.667 

.333 

0 

Mounted  particles. 

1 16(a) 

29 

Zr-U 

.667 

.333 

0 

Mounted  particles. 

116(b) 

55 

Zr-U 

.705 

0 

0 

Mounted  particles. 

117(a) 

64 

Zr-U 

1.000 

0 

0 

Mounted  particles. 

117(b) 

65 

Zr-U 

.667 

0 

.333 

Mounted  particles. 

117(c) 

66 

Zr-U 

.667 

.333 

0 

Mounted  particles. 

117(d) 

HI 

ZrH-U 

4.400* 

0 

0 

Mounted  particles. 

117(e) 

H2 

ZrH-U 

4.384* 

0.016* 

0 

Mounted  particles. 

117(f) 

H3 

ZrH-U 

9.800* 

4.900* 

0 

Mounted  particles. 

Pressures 

measured  in 

psia  instead  of  atms. 
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APPENDIX  II 


ANALYSIS  OF  CHEMICAL  REACTIONS  OF  A  HYDRIDED 
FUEL  ELEMENT  DURING  RE-ENTRY 

i 

An  analysis  of  the  behavior  of  a  zirconium -uranium  alloy  in  a  high  temperature,  air 
atmosphere  environment  shows  that  it  oxidizes  according  to  the  parabolic  rate  law  and  absorbs 

I 

both  nitrogen  and  oxygen  into  it.*9'34'*  The  oxide  formed,  which  is  ZrOa,  also  absorbs  oxy- 

| 

gen  and  nitrogen,  as  is  evidenced  by  the  phase  diagrams.49  Initially,  it  is  assumed  that  small 
alloy  additions  of  uranium  (say  up  to  4  atomic  percent)  do  not  change  the  fundamental  be¬ 
havior  of  zirconium  in  an  oxidizing  atmosphere.  The  data  will  be  shown  to  be  in  essential 
agreement  with  this  assumption. 

A.  FACTUAL  INFORMATION 

The  data  available  for  analysis  consist  of  the  following  known  items. 

1  .  Weight  gain  measurement*10 continuous  in  time,  on  sample  fuel  element 
material  specimens  for: 

a .  100%,  50%,  and  10%  air  atmosphere  environment  at  one  atmosphere 

I 

of  pressure .  The  dearth  of  air  pressure  in  the  10%  and  50%  air  cases 
was  made  up  by  the  addition  of  argon  gas.  The  argon  is  known  not  to 
react  with  the  fuel  element  material . 

b.  500°C,  800°C,  1100°C,  and  1400°C  temperature  environments  were 
used  with  each  of  the  atmospheres  listed  in  a. 

2.  The  starting  compounds  and  reaction  product  for  each  of  the  weight  gain 
experiments  are  known.  These  include: 
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a.  Initial  composition  or  weight  percent  of  hydrogen,  zirconium, 
and  uranium. 

b.  The  final  composition,  hydrogen,  nitrogen,  oxygen,  zirconium, 
and  uranium. 

c.  The  coating  formed  on  the  sample  specimen  surfaces,  i  .e.  ZrOa  . 

3.  The  length  of  time  required  for  the  fuel  elements  to  reach  an  essentially 
uniform  temperature  under  the  environmental  conditions  described  in  lb 
above  is  known. 

4.  The  rate  of  hydrogen  evolution  from  the  Zr-H-U  material  as  a  function  of 
the  temperature  for  the  range  given  in  lb  above,  when  heated  in  an 
argon  atmosphere,  is  known. 

5.  The  quantity  of  hydrogen  that  the  fuel  element  material  will  maintain  in 
an  activated  absorption  state  is  essentially  independent  of  the  external 
hydrogen  pressure,41  varying  roughly  from  the  one-eighth  to  the  one- 
fourteenth  power  of  the  external  pressure  at  a  given  temperature. 

6.  The  rate  of  oxidation  is  essentially  independent  of  the  external  oxygen 
pressure,  so  long  as  oxygen  is  supplied  at  the  outer  surface  of  the  oxide- 
nitride  layer  (which  is  formed  on  the  material),  as  fast  as  it  is  used  up  at 
the  interface  between  this  layer  and  the  fuel  element  proper.37 

7.  Phase  diagrams  for  the  activated  absorption  of  oxygen  in  zirconium  and 
zirconia  are  available.46 

8.  Insofar  as  the  heat  energy  released  due  to  the  chemical  reactions  of  oxygen 
and  nitrogen  with  the  fuel  element  material  is  concerned,  the  energy  pro¬ 
duced  per  unit  weight  gain  of  the  fuel  element  material  is  essentially  the 
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2ZrN. 


same  for  the  two  reactions,  Zr  +  ^  — ♦ZrO*  and  2Zr  +  Nt  — > 

These  are  the  principal  reactions  to  be  considered.37  An  X-roy  analysis 
of  the  surface  material  produced  on  the  fuel  element  material  under  simu¬ 
lated  re-entry  conditions31  and  during  the  weight  gain  experiments10  shows 
that  the  surface  material  is  ZrOa  . 

9.  Zirconium  and  dilute  binary  alloys  of  Zr  with  other  metals  will  oxidize  much 
more  rapidly  in  an  atmosphere  containing  a  few  percent  water  vapor  than  in 
a  dry  atmosphere  .34 

10.  When  the  fuel  element  material  is  exposed  to  a  heating  environment  similar 
to  that  that  it  would  experience  during  re-entry,  but  containing  an  excess  of 
water  vapor,47'48  the  sample  fuel  elements  being  tested  partially  me  I  red  and 
ablated  away  during  the  tests. 

1 1 .  Samples  of  fuel  element  materials,  when  subjected  to  varying  arc-image 
heating  rates  and  atmospheric  pressures  that  simulate  those  that  would  be 
experienced  by  a  fuel  element  along  a  variety  of  re-entry  trajectories,  do 
not  lose  their  geometrical  integrity  .31  However,  they  do  lose  hydrogen  and 
do  have  coatings  of  Zr02  formed  on  their  surfaces. 

12.  The  solubilities  of  oxygen  in  Zr  and  ZrQj  are  approximately  equal* 

13.  When  heated  in  an  air  atmosphere  for  long  periods  of  time,  the  fuel  element 
material  exhibits  "break-away"  oxidation.  Due  to  the  large  thicknesses  of 
the  oxide  involved,  this  is  not  attributed  to  a  mismatch  in  thermal  expansion 
coefficients  between  the  alloy  and  the  oxide. 

14.  It  can  be  shov/n  analytically  that  the  most  rapid  possible  rate  of  oxidation 

(or  other  chemical  reaction)  in  an  ion  diffusion  controlled  process  is  thatrate 

described  by  the  parabolic  oxidation  law. 
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B.  MATHEMATICAL  MODEL  FOR  CHEMICAL  REACTIONS 
1 .  Correction  Factors  for  Weight  Gain  Measurements 

The  weight  gain  measurements  made  by  ARF,  described  in  Item  1  above/  were  made 
on  short  cylindrical  samples  (0.585  cm  long,  1 .905  cm  diameter).  Because  the  oxide  layers 
formed  on  the  surface  of  these  test  samples  were  thick,  and  because  the  oxidation  rate  laws 
are  based  on  weight  gains  per  unit  area,  it  is  necessary  to  correct  the  weight  gain  data  for 
the  fact  that  the  surface  area  of  the  unreacted  material  decreases  as  the  weight  (thickness  of 
the  oxide  layer)  of  the  sample  increases.  By  making  use  of  the  geometry  of  a  cylinder,  it  can 
be  shown  that  the  area  of  the  sample,  A(t),  that  remains  at  the  time,  t,  and  hence  is  available 
for  oxidation  is  simply  given  by 

A(t)  =  2n  (R  -S)  [(R  -S)  +  2(L-S)]  (91) 

where  R  =  initial  cylinder  radius, 

L  =  one-half  of  the  initial  length  of  the  cylinder,  and 

S  =  thickness  of  the  metal  that  has  been  converted  to  a  coating  (i  .e. ,  an 
oxide  and/or  nitride). 

It  can  also  be  shown  in  an  analytical  manner  that  if  f(t)  is  the  fractional  weight  gain  of  the 
sample  at  the  time,  t,  and  K  is  the  weight  of  the  reacted  compound  containing  one  gram  of 
the  metal,  then  the  value  of  S  is  the  solution  to  the  following  relation: 

R3  L  fftl 

S*  -  (L  +  2R)  S3  +  (R3  +  2RL)S  -  p  =  0.  (92) 

The  total  weight  of  the  sample  at  time,  t,  is  the  initial  weight,  We,  plus  the  weight  WQf(t). 
The  data  on  weight  gain  measurements  were  taken  on  a  strip-chart  which  showed  the  integral 
weight  gain  as  a  function  of  time.  Because  the  relative  magnitudes  of  the  possible  reactions 
that  took  place  during  the  data-gatherlng  time  interval,  these  data  must  be  reduced  to  finite 
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difference  form  for  the  proper  interpretation  of  the  results.  The  weight  gain  per  unit  area 

during  the  i_th  time  interval.  At;,  can  readily  be  shown  to  be  given  by 
,A  .  _  Af(At;)W0 

Am,  (At;)  - -  (93) 

A(t') 


where  Af(At;)  is  the  change  in  the  fractional  weight  of  the  sample  during  the  time  interval,  Atj , 

and  t1  =  t  At; .  It  is  clear  that 
2 

I  Atj<  =  t,  (94a) 

k=  1 


Z  Am(Atk)  =  Anis(t),  and  (94b) 

k  =  1 


Z  Af(A  tk)  =  f(t) .  (94c) 

k=  1 

These  equations  are  applicable  to  any  cylinder  that  is  undergoing  a  weight  gain  due  to  a  surface 
reaction  that  forms  a  coating  of  products  on  the  surface  thereof.  The  quantity  f(t)  is  interpreted 
as  the  weight  gain  (or  loss)  due  to  the  formation  (or  decomposition) of  the  surface  coating. 


2.  Description  of  Reactions 

The  weight  gain  (or  loss)  of  the  element  is  due  to  oxidation  and/or  nitridation  (Am,), 
to  nitrogen  and  oxygen  absorption  (Ama),  and  to  dehydriding  (Am^),  so  that  the  total  weight 
change  at  a  given  time  is 

A  m(t)  =  Ams  +  AmQ-A  mH  •  (95) 

The  weight  gain  data,  when  corrected  in  accordance  with  eqs.  (93)  and  (94),  lead  one  to  sus¬ 
pect  that  at  first  there  is  some  loss  of  hydrogen  and  absorption  of  nitrogen  and  oxygen,  together 
with  oxidation.  After  some  time  has  lapsed,  the  process  follows  a  true  parabolic  law,  and, 
finally,  for  the  very  thick  coatings,  exhibits  breakaway  oxidation.  In  view  of  the  fact  thot 


oxygen  and  nitrogen  have  almost  identical  atomic  sizes  and  masses,  no  distinction  will  be 
made  between  them  insofar  as  mass  increases  due  to  their  absorption  into  the  sample  material 
during  weight  measurements  are  concerned.  With  this  simplification,  and  in  view  of  Item  8 
above,  eq.  (95)  can  be  written  as 

Am(T,  t,  p)  =  [Ast]  8  e  Q,/2T+  Aq  p^8te  E^T-  Amn(T,  t).  (96) 

Here,  RgQ  and  RgE  are  activation  energies,  when  Rg  is  the  universal  gas  constant,  A$  and 
Aq  are  constants,  and  p  is  the  partial  pressure  of  the  component  of  the  gas  being  absorbed. 

The  first  right-hand  side  term  in  this  relation  represents  the  parabolic  chemical  reaction  law 
and  the  second  the  activated  absorption  law.36 

a .  Determination  of  activated  absorption  constants.  —  The  values  of  Aa  and  t  are 
determined  as  follows.  Holding  T  and  t  constant  and  letting  p  vary  for  two  sets  of  t  and  T, 
eq .  (96)  shows  that 

Am(Tx  ,  tx  ,  px)  -  Am(Tx ,  tx,  pa)  =  Aatx  e  1  [px^2  -  p/*]  (97a) 

and 

Am(Ta,  t8,  px)  -Am(Ta,  ta,  p,)  =  AQta  e"^Ts£pl*»-  p2^a]-  (97b) 

These  two  equations  determine  the  values  of  Aq  and  E.  The  values  of  t  to  be  used  must  be  suf¬ 
ficiently  small  to  assure  that  the  absorption  process  is  not  seriously  inhibited  by  the  oxide  coat¬ 
ing  that  is  being  built  up  on  the  surface  of  the  sample.  Examination  of  the  weight  gain  data 
show  that  Ama  is  only  a  few  percent  of  Am$,  being  less  than  1 .7%  weight  gain  in  all  cases. 

Its  values  are  essentially  masked  by  errors  in  the  data  and  hence  can  be  considered  as  negligible 
insofar  as  the  present  analysis  is  concerned.  The  oxygen  and  nitrogen  absorption  effects  would 
be  even  less  under  actual  re-entry  conditions  because  of  the  low  pressure  in  the  boundary  layer 
about  the  fuel  element  during  the  initial  stages  of  oxidation. 
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b.  Determination  of  hydrogen  loss  rote.— By  heating  the  fuel  element  material  in  an 
argon  atmosphere  and  noting  the  rate  of  hydrogen  loss,  it  is  immediately  apparent  that  the  loss 
rate  is  an  exponentially  decreasing  one.  This  suggests  the  following  differential  equation  : 

d  M^r 

“dt — =  (MHo-MHrH, 

where  Mj^r  is  the  amount  of  hydrogen  that  has  been  released  per  unit  area  from  the  sample  sped- 
ment,  MHo 's  the  anrount  of  hydrogen  that  was  initially  available  for  release  (per  unit  area)  at 
the  temperature  in  question,  and  X  is  a  temperature  dependent  parameter.  Integrating  this  ex¬ 
pression  and  evaluating  the  constant  of  integration  by  making  use  of  the  boundar’  condition 
Mnr  =  0  at  t  =  0  gives 

M|^r  =  M|-|0  0  “  ®  )  *  (98) 

Making  use  of  the  dehydriding  rates  reported  in  reference  10,  a  plot  of  X  versus  T  can  be 
obtained.  This  plot  is  a  straight  line  whose  equation  is 

X  =  1.01  x  lO^T  -  1.92.  (99) 

This  relationship  is  valid  for  T  >  1900°R .  Below  this  temperature,  the  loss  rate  can  be  considered 
to  be  zero  for  the  practical  problem  at  hand . 

The  hydrogen  that  is  released  diffuses  through  the  oxide  layer  that  is  being  formed  on  the 
surface  of  the  fuel  element  material.  The  rate  of  loss  through  the  oxide  layer  per  unit  area  is 
governed  by  a  diffusion  equation  of  the  form 


d% 

dt 


-  a  grad  C^, 


(100) 


where  is  the  amount  of  hydrogen  that 


has  diffused  through  the  oxide  layer  and  escaped  into 
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the  surrounding  environment,  C j_j  is  the  hydrogen  concentration,  and  a  is  the  diffusion  con¬ 
stant  multiplied  by  thedensity  of  the  oxide  film.  Since  the  concentration  just  outside  the  oxide 
surface  is  essentially  zero,  it  follows  that 

grad  CH  =  Mhr  I  Mh^  .  (101) 

Amj 

Substituting  eqs.  (98)  and  (101)  into  (100)  gives,  upon  rearrangement. 


d  Mh-1  _  a  MHt  _  a  MHoO  ~  e~Xt) 
dt  /T^T  j  kst  ' 

where  ks  is  the  parabolic  rate  constant,  i.e. 


002) 


-Q/T 

ks  =  As  e  .  (103) 

This,  then,  is  the  differential  equation  describing  the  rate  of  hydrogen  loss  (per  unit  area)  through 
the  surface  of  the  oxide  coating. 

To  solve  eq.  (102),  note  that  it  is  of  the  form 


-/■  +  f(x)y  =  g(x) 
ax 


(104a) 


The  complimentary  solution  (solution  to  homogeneous  equation)  is  simply 

_  r  -JfM  dx  r  _ 

y  =  C  e  ,  C  - 


constant. 


(104b) 


letting  C  =  C(x),  and  using  the  method  of  variation  of  parameters,  gives,  upon  differentia¬ 
tion  and  substituting  into  (104a), 

Jf(x)d 


CM  =  J[gM  + A 

where  A  is  the  constant  of  integration.  Substituting  eq.(104c)  into  (104b)  gives 


(104c) 


y(x)  =  exp 


[-  J  f(x)  dx  ]  ’  A  +  J  g(x)  exp[  J f(x)  dx  ]  dx  j.  . 


(104d) 
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Substituting  eq .  (102)  into  this  general  solution  gives 


-h  =  mh*H [Ji-rOH 


005) 


Now,  %  (a/-/ kg  t')  dt  is  readily  integrated  to  give  2a(t/ks)%,  thus,  (105)  becomes 


_  --  2a  /t/ks'  a  M|Ho  2a  /t/ks' 
AmH  “  Ae  -  7= — ~  e 

V  *s 


I 


-(X  t  +  2a  J  t/kg  ) 


-  2a  /  t/ks  ^ 


■dt 


(106) 


The  las*  integral  in  eq.  (106)  is  readily  integrated,  that  is, 

»  e-2a  v/t/ksV  -2a  J t/ks ' 

— - — ; - dt  =  - —  e 

J  /  t  a 


If  q  =  t^*  is  substituted  into  the  first  integral  in  eq.  (106),  it  becomes 

-(\t  +  2a/~t/ks) 

- - dt  =  2 

/r 

Then,  let  u  =  Xqa+  2aa/J~V^ ,  to  obtain 


e-(Xqs+2aq//i5)dq. 


2je'°  _ - jg 


Integration  by  parts  gives 


(107a) 


(107b) 


Now,  change  variables  by  letting  v8  - 

- 2 1  du 
Integrating  this  expression  by  parts  gives 


in  the  integral  of  eq.  (107b) 
Je'Xvl(2vdv). 


to  obtain 
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The  bracketed  expression  is  simply  erf^v/x"*).  Substituting  this  into  (107b)  and  changing 
variables  back  to  t  gives,  upon  simplification, 

2(/8+  a/x/kT )  e"(Xt  +  2a  ^t/k»  >  [l  +  e“2aa/Xks]  -  2  e"a3/Xkserf(Xt  +  a//TkJ  )  (107c) 
Substituting  eqs .  (107a)  and  (107c)  into  eq.  (106)  gives  the  complete  primitive  of  eq.(102)  as 


AmH 


=  _e2a/i7kT+ 


_  I 

-2a2Aks]  -  2  /2tT  e"0*^5  erf(/Tt*  +  o//xkT  )  +  /^L  e”2a  7  (108) 

•  a  j 

The  constant  A  can  be  evaluated  by  making  use  of  the  boundary  condition  that  ftm^  =  0  at 
t  =  0.  This  gives 


A  = 


MHc 


7k  s 


(2a/\/ET) 0 +  e2a3/Xks)  - 2 /2tT  e’^/Xtcs  erf(a/XkT 


+  /k^a)^ 


(109) 


Substitution  of  this  value  of  A  into  eq  .  (108)  gives  the  amount  of  hydrogen  lost  per  unit  area 
in  a  time  t  as 


,  x  _  qMh°  2a/t7kT 
6mn(0  _  ® 


0  +  .'2aVXk«)  +  2  iw  (AT 


+  q/J\ kj’)  -  *rf(a//\k|’  J  +^~  (1  -  •  2a ' ) 

\ 

4  2(rr  4  oa/ct  >  0  4  o-2a*Ak>)  .‘(u+ 2a/^7>  ■ 


(110) 
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The  diffusion  parameter  must  be  determined  experimentally.  This  is  accomplished  by 
measuring  the  quantity  of  hydrogen  left  in  the  sample  specimens  after  they  have  been  exposed 
to  the  weight  gain  tests.  There  was  quite  a  spread  in  the  initial  weight  percentages,  Mhq,  of 
the  test  sample  specimens,  varying  from  approximately  1  .4  weight  percent  to  1  .7  weight  per¬ 
cent.  Because  of  this  nonuniformity  among  the  various  specimens  tested,  the  mean  value  of 
MHo  =  1  .55  weight  percent  was  taken  for  determining  a.  Various  sample  specimens  that  had 
been  exposed  to  the  weight  gain  tests  were  analyzed  for  hydrogen  content  by  vacuum  fusion 
analyses.  This  results  in  a  corresponding  spread  in  the  values  of  a.  The  value  of  a  so  deter¬ 
mined  is 

a  =  Ah  e"  Qh/T,  where  Ah  =  0. 10,  and  Qh  =  1 -B4  x  104  .  (Ill) 

Substituting  these  values  and  eq.  (99)  into  eq.  (102)  gives  the  rate  of  hydrogen  loss  per  unit  area 

^[AmH(0]s-^2tM--  -aMHo[l  -  exp  (1.92  -  1.01  x  10JT)]  .  (112) 

The  value  of  AmH(t)  is  given  by  eq.  (112),  and  the  value  of  MHo  has  to  be  known.  For  the  par¬ 
ticular  samples  used  in  the  weight  gain  tests, 

W 

MHo  =  1  *55  x  10-  — —  =  1 ,69  x  10“3grams/cm8  . 

A(o) 

The  values  of  A$  and  Q$  are  to  be  taken  from  eq.  (116)  below. 

c.  Determination  of  the  parabolic  rate  constant,— When  the  weight  gain  data,  as 
modified  by  eq.  (93),  are  plotted  on  log-log  paper,  the  result  is  a  family  of  parallel,  straight 
lines  for  t  greater  than  about  15  to  20  minutes,  one  line  for  each  temperature  (see  Fig .  7)  It 
is  readily  determined  that  the  process  involved  is  parabolic  in  nature,  that  is,  each  of  the  lines 
car.  easily  be  fitted  by  an  equation  of  the  form 


86 


Amf  (T,  t)  =  k(T)t  +  C(t) . 


013) 


By  holding  T  constant  and  varying  t,  one  can  calculate  the  value  of  k(T),  i.e. 

A  mfCTjtj)  -  Anrif (T ,  t8)  =  k(T)  (tx  -  t8).  (114) 

Here,  it  has  been  assumed  that  C(t),  which  arises  from  Amg  and  Am^  is  actually  constant  for 
large  values  of  t,  say  t>20  minutes.  The  fact  that  the  family  of  lines  are  truly  parabolic 
supports  this  position.  That  is,  if  there  are  changes  in  Ama  and  Am^  for  large  values  of  t,  they 
are  overshadowed  by  Ams  to  within  the  limits  of  the  accuracy  of  the  data . 

The  reaction  rate  constant  is  known  to  obey  Arrhenius'  equation,23  that  is, 

k(T)  =  Ae'Q/T.  (115) 

The  value  of  k  is  determined  from  eq.  (114),  using  the  temperatures  1390^11,  1930^R,  2470°R, 
and  3010°R  .  There  is  a  phase  change  in  the  fuel  element  material  at  approximately  2040PR. 

This  phase  change  changes  the  value  of  k.  If  log  k  is  plotted  versus  1/T  for  the  1390°R  and 
1930°R  points  and  for  the  2470°R  and  1930°R  points,  the  pair  of  straight  lines  connecting  these 
two  pairs  of  points  cross  at  the  phase  change  temperature  (see  Fig.  8).  This  is  the  result  ex¬ 
pected,  and  it  lends  even  greater  confidence  to  the  values  of  A  and  Q  so  determined.  They  are 

As  =  0.356  (gms2/cm*-min),  Qs  =2.59x10*,  for  T  ^  2040°R ,  (116a) 

As  =  92.5  (gmsJ/cm*-min),  Qs  =  3.79x10*,  forT^  2040°R.  (116b) 

With  these  parameters  determined,  the  mass  increase  at  any  time  can  now  be  determined  by 
eq.  (95)  for  any  combination  of  the  sets  of  variables  (T,  t).  Also,  the  rate  of  hydrogen  evolution 
is  determined  by  using  the  results  of  eq.  (116)  in  eq.  (112). 
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